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l!he launch of I£iNDS2^ has for the first tints provided the water rer 
source manager and practical hydrologist with broad proai«ots for 
efficient acgoisiticn of essentially real**time data. Ihese are usable 
for hydrologic land use assessment, surface cover classification, phy- 
siographic analysis, surface water inventory, and for the retraction 
of information pertinarit to soil properties. This infemation has 
value rot only by and in itself, but also to construct the watershed 
transfer function for i^^drologic planning itodels aimed at estimating 
peak outflow from rainfall inputs. 

The reduction of satellite data to practical, operational information 
requires a clear, easily applicable methodology for oaiverting these 
data into quantitative hydrologic pairameters. 

The fundamental objective of this effort is the development of such a 
methodology and its transfer to hydrologic users. It was realized 
that such technology transfer could be made far more effective ky the 
parallel development and eventual demonstration of the results of a 
model, specifically structured to t^e full advantage of the capabil- 
ity of UMDSM — for example, its frequent recurrence and cemsegusnt 
ability to detennine seasonal variations in the watershed's conditions. 
The category of Planning Models vas diosen for development and demcn- 
stratioji because of its great practical irtport^ce in the design of 
waterworks, because of the wide diffusion of such models devm to capil- 
lary levels within tl^ hierarchy of water resources users, and because 
their inplamentation is relatively siitpler than real-time management 
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inod^Sf thus making optimum use of the resources available for this 
effort, 

Oonseguently, the effort wts structured aloi^ ti-ra major routes: the 
defoeloipment of a hydrologic planning model specifically based i:pan 
renotely sensed ii^uts, including its test and verification from ex- 
isting records; r.id the applicc^tlon of IJ^IDSM' data to st^plying the 
mode's quantitative parameters and ooeffipients. Included was the 
investigation of the use of XJM:)SKr data as information inputs to 
all categories of hydrologic models requiring quantitative surface 
parameters for their effective functionij^, 

"®ie effort idius far has consisted of two phs^ea, Ihe first focused 
on the d^initlon of the "drivers" - those hydrologic processes to 
which peak runoff is most sensitive - arui upon the synthesis of a 
sisiple yet effec^ve model for the estimation of Icaig-recurrence out- 
flows, Ihe results of the first phase effort were presented in the 
Fii'al Kfiport, "Ihe JJ^lication of Ranote Sensing to the Develcpnent 
and Ponnulation of ideologic Planning Models," dated January, 1975. 
Ihe seocoid phase has extended this work to include the development 
of a routing model for use in sensitivity analyses, and a quantita- 
tive investigation of the accuracy and ocmpleteness of the hydrolo- 
gic information which can be extracted fresn remotely-sensed imagery. 

This document reports the findings and oonclimions of the Phase-two 
effort: it includes a summary of the results of the earlier work. 




OffiKEER 1 


mjmi OF THE FIRST EEFOBT 

Of critical cxanosm to vrater resources planners and engineers is the 
ability to forecast ]^ak flcftv events. The capacity to estimate the 
magnitude and duration of largerrecurrenoe cFotflows has a significant 
in^ct i^xsn the accuracy of sizing and designir^ watervrorks, and thus 
oaa their cost. 

The tool available to the planner for these purposes is the hydrologic 
model, Although the ii^nita of different models vary/ all re^iire sig~ 
nificant quantities of p^sjogrc^phlc and hydrologic InfomBtion; these 
data axe typically e}^>ensiv3 to obtain and are ofte;) only paartially 
atvailable, Remote sensing offers a as» source of information whicdi 
foxm^ly had to be acquired by less efficient means or ignored al- 
together. 

The first phase of this effort conducted fron February to December 
1974 addre^ed four pertiTient topics : 

1) Ideratificeition of the "drivers" of peak flow events, i.e., 
the l^dxologic phencmena (infiltration, antecedent soil 
moisture, etc,} to which the watershed’s outflow is most 
sensitive. 

2} The dsvelopnent of a model oompatible to the maximun degree 
with renotely-sensed hydrologic ii^uts* 

3} Verification of the model for actual watersheds. 


4) Pxeliiniiiary idsRtification of the efficiency of remote sen- 
sing in detemining the parameters of the model. 

X.l) Innnsstigaticm of l^river Phenooena 

Ihe purpose of this investigation was to achieve a valid statistical 
OGB^parison of the rates and magnitudes of the hydrologic processes 
oontributing to the runoff fron long-recutrenoe events. OMs was ac~ 
ccnplii^ed over a significant sanple of watersheds , with wide var- 
iaticm of cdimatology, terrain, and physiograpJ^'. Ihe ccai5>arison al- 
lowed the deterndnation of whidi are important and which can he neg- 
lected without significant loss of a<xmracy» 

Bain falling on a waters!^ is subject to several processes vMch ab- 
stract water and govern the flow* IChose which produce most sig- 
nificant changes to flow rates and volumes "drite*' the basin outflow. 
Table 1 describes the i^drologic processes; Table 2 presents the 
‘•drivers** of each. Figure 1 synopsizes their relative contributions 
to the peak event. It shows that several processes can be omitted in 
the fbnRUlation of a p^c rate model because of their limited impact. 
For eKample, the rates oorrespcmiding to interflow, percolation, and 
ev^otranspiration are very slow in ocsiparison with other processes 
such as rainfall, infiltration and overland flow. Also, interception 
and depression storage beccme saturated early in large rainfall events. 

Therefore, except for very special circumstances, peak flew can be 
adequately modelled by considering only precipitation, infiltration arid 
surface flow - both overland and in the channels. 

)Y THB 
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5U?f!lARy DESCRIFnONS OF HZDROIiOGIC PB0CESSE3 


HXDROLOGIC PFSJCESS 

DESCRIPTION 

Interception 

Moisture caught and' stored on plant leaves and 
steans or other inpermeahle objects; eventually 
j evaporate back into the atitiosphe3?e.' 

]jnflltration 

Downward movanent of water from the surface into 
the soil. ■ 

: .. V „ 

A) liitepflow 

f. Lateral subsurface water movement towai^ stream 
channels , 

B) Percolation 

Downward movement of water thoKsu^ soil to 
groundwater (area where pores of soil or rock 
are filled- with water) , 

C)' Base Kunoff 

Water from Interflow and percolation which moves 
underground to the channel. 

Evapotranspiration 

Upward movement of water in ^seous state from 
the surface. 

A) Evaporation 

B) Ti?anspiration 

, Movement of water tnrough plants to the 
atmosphere, ■ • 

Precipitation Excess 

Retention of excess rainfall in surface depres- 
sions, 

A) Depression 
Star*age 

- B) Surface Plow 

Uninfiltrated water which flows' over land , 

surface to stream channels, . * 

C) Channel Plow 

Plow of water in natural channels, . 

Total Runoff 

Sum of runoff from underground -processes (base 
runoff) and overland flow (direct runoff) . 



TABLE 2 


POTENTIALLY IMPORFAm DHIVERS AS RELATED TO BKDBOLOGIC PROCESSES 


H2DRDL0GIC PROCESS 

PRINCIPAL DRIVERS 

SECONDARY DRIVERS 

Overland Flow 

\ ■ 

1 I III 1 1 1 1 .1 1 

Slope 

Roughness of Soil 
a Cover 

Drainage Density & 
Pattern 

— — 

Infiltration 

Soil Pertneability 
Antecedent Soil 
Moisture 

Soil Moisture Capa- 
city 

Vegetative Cover 
Slope 

Water Turbidity 
Tenperature 

A) Int^low 

Soil Pemeabllity 
Subsurface- Itolsture 
Gradient 

Plow Length, Slope 

4 ^ 

B) Percolation 

Soil Pemeabllity 
Subsurface jfoisture 
Gradient 
• Soil Depth 

— 

Evapotransplration 

Tenperature ’ 
Antecedent Soil 
Moisture 

Soil Permeability 

Water Turbidity 
Wind 

A) Evaporatlcai 

B) Trar^lratlon 1 

Tsrperature 
Solar FMiation 
Vegetative Cover 
Antecedent Soil 
Moisture 

Wind 

Depression Storage & 
Detention 

Depression Density 
Cover Retention 

Slope 

Inta?ception 

i 

Duration of Itelnfall 
Intensity of Rain- 
fall 

Cover Coirposition, 
Age, Der^lty 

Evaporation Rate 
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^Ehe set o£ the principal drivers of planning models can thus be de* 
fined as shown in Table 3 which also indicates that several of the 
drivers are raiotaly measurable or inferrable. The means by which 
this can be acoon^lished be discussed in Chapter XV. 

1,2) IDev^opnent of Bemote Sensing model 

The £ollo\<dng cnriterla were foUot^d in the structuring of a peak-rate 
nodelt 

1] The model would be modular, to allow the user flexibility of 
application. 

2) The model wuld provide the long-recurrence peak outflow rate; 
development of a model to yield the hydrograph was reserved 
for a later phase of the effiort. 

Sensitivity analyses of the sub-surface abstraction module were per- 
formed on an ERI 380 Hybrid oaiputer to detamine the effects on run- 
off of t^e. sub-surface hydrologic drivers ; soil itoisture content and 
capacity, and soil perm^bility. Figures 2 and 3 exen^ify two of 
the cnnpiter runs. 


^Ihe pe^ rate model was constructed from t!^ following nodules: 


Rain Recnjrrence Module; 


an enppirlcal formulation of the form 
(t + d)“3 


( 1 ) 


i = rain rate, nv'sec, 

T t= recurrence period, years 




. . . ■■ OF ‘FF)'- 

paSe is poor 


where: 
























Rote,cm/*w' 



I01L UPON RUNOFF 


^ LOM., H»«* 

COUNTY, MD. 


upper Layer Depth- 03m 
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t S' rain duration, hrs, 
a 2 > a^f <3 = constants, functions of location 

This was used to calculate the expected magnitude of large re- 
currence rainfalls. 

Rainfall Spatial Correction Module ; the introduction of a factor 
to convert point rainfall to its areal equivalent for large basins. 
As watershed area increases^ effective rainfall must be diminished. 
An emirical relation, as shown in Figure 4, was used for this pur- 
pose. 


Sub-surface Abstractions (Infiltration) Jfodule; the Holtan inr 


filtration ©guation; 

i = GI • a • (S^ - I)^'^ + Ig (2) 


m 

where: I 

GI 


a 


I 


I 

f 


=5 infiltration rate 

= maturity of cover index 

= average vegetative cover factor 

= average available soil moisture 
capacity 

= cumui.ative infiltration 
= saturated infiltration rate 


Holtan* s equation was selected for use because 1) it e^licitly in- 
cludes surface observables (the a and (31 factors) which are poten- 
tially rennotely sensible, 2) it is currently being applied over a 
diverse range of cover and soil t^^pes so o:tipirical verification ex- 
ists and 3) its results do not differ appreciably frcati those of 
other widely-used fomulations. 




Overland Flow Module ; a formulation viiich relates maximum over- 
land flo»; rate to watershed slope, area, surface friction, channel 
length and rainfall input. 


vdiere: 


Q * 


2U C 


3 

1 

<fi (D (n) 5 

0 4 — 

2 3 

* 

e 5 S 10 (3600) 



(3) 


3 2 

Q = maximum out flew rate, m /sec/lon 
L channel length, m 
1 = length of average strip, m 
n"= average Manning's "n" 

(f> - routing factor 
5 = K T “2 


K = infiltration and spatial oorreettion factor 
T = rain recurrence interval, years 
s - average slope, nv^ 

02 , = constants, function of location 


1.3) Verification of the Model 


Validation of the model required testing against a sat of real watersheds 
possessing long-term records and r^3<resenting a variety of environmental 
conditions. Such a set has been developed by the U.S. Department of 
Agriculture, Agricultural Research Service (ARS) . Prom this set of 
^pcxndmately 250 watersheds, 158 basins with area greater than 40 hec- 
tares vere' selected to form the analytic sanple. Kis fifty-year re- 
currence outflow was selected as the test event. Initially, the fifty- 
year event was chosen because it represented a realistic estimate of 
water resources project life and tecause it was consistent with the — 
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lengbh of available records. The selection vras made sxibject to further 
verification which has been oonpleted and is described in Fection 3.3. 
Initial verification of the model was nade on a s st of nine basins 
selected for geogre 5 )hic and l^rologic diversity. A typical descrip- 
tion of such a test basin ^ip^rs in Figure 5. The parameters necess- 
ary for application of tbe model were calculated fron the AES data, A 
stminaxy of the parameters for the test subset of nine basins is pro- 
vided in Table 4, 

The coiputed flow rate and the rate statistically derived fron the 
measurement records were ocnpared. Oonputation of peak flew was also 
mads using other models in oantron usage - the Eational Formula, the 
S.C.S, model, and Oodk's equation - for the same watershed san^le. 

The results are presaited in Table 5 and Figure 6. The peak-rate 
model yielded estimates within 1 . 15% of those derived fron runoff 
records for seven of the nine test watersheds. 

T!he results for two watersheds - Raynol<fe, Idaho, and Chickasha, 

CJklahatB, exceeded the ±15% accuracy bound. It ^>peared that these 
large errors might be attributable to the "ccRplexity” of the two 
basins, both of which are conposed of numerous sub-fcasins of diverse 
oharacteristics^henc^ requiring more ccaxplex routing than inoorpora- ■“ 
ted in the model. The development of the routing module will be de- 
scrilaed in Chapter III. 

1.4) Identification of the Eole of Eemote Sensing in I^drologic 
Modelling 

A visual analysis was performed of one watershed at ChieJeasha, Oklahoma, 
fron ii^ack and white. Band 5, LAKESAT imagery taken on October 20, 1973, 


FI6URE 9 


RIESELtWACO) TEXAS ECO- 7 


Area 4.5 Km^ 

Slope .02! m/m 2 

Length of Channel 6.6 Km 

Drainoge Density 1/680 m/m' 
(t+i) -78 


Cover 

60% Pasture 
6% Small groin 
3% Corn 
7% Cotton 
9% Row grains 
2% Gravel 6 paved roads 
1 3% Other, mostly weeds 










Soils 

66% Wilson cloy loam 
24% Burleson Heiden clay 
4% Frio cloy loom 
3% Crockett loam 
2% 6urles(m cioy 
I % Houston Slock clay 






REPRODUClBlLm^ OF ' 
ORIGINAL PAGE IS P('\ 



TABLE 4 

PEAK RATE MODEL PARAMETERS 


LOCATION 

fceo 

MAIN CKW' 
NEL LEN6TH 
L,m 

OVERLAND 
FLOW LEN81H 
1.0- 

(D 

3 

i 

*w 

k3/H^IW0 

K 

OANVItXE,VT. 

4.»«ltf* 

9314 

808 

2.0 

.079 

.12 

.91 

.28 

COSHOCTON, O. 

l.4*Kr® 

518 

163 

1.73 

.048 

.172 

28.5 

.71 

BLACK8BUNB,VA. 

1.4* Kf® 

8008 

167 

12.1 

.039 

.129 

1.01 

.62 

OXFORD, MISS. 

i.4*icr® 

12,000 

898 

2l8I 

.060 

.114 

iOlS 

.6) 

FOmiMORE, Wise. 

14*10*8 

1849 

166 

926 

.084 

.080 

12.8 

.70 

CHiCXA8KA,OKLA. 

2J0*©*® 

94719 

800 

94.4 

.038 

.088 

.08 

j 

.68 

WACO TEX. 

2.2*10*® 

3288 

283 

3.1 

! 

.038 

.021 

lt.9 

.76 

SAFFORD,ARIZ- 

... ...j 

i.2* lor*® 

8549 

133 

6.6 

.020 

j 

.020 

53 

.73 

i 

RSYHSLDS 9 QAHO 

j 

2,8* 10*^ 

21481 

189 

! 

32.2 
! 

.096 

.176 

.001 

.60 










TABLE 5 


FC 

FL 


1. OonviHe Vt. 

2. CoshoGtontOhb 

3. Blacksburg -Va 

4. Oxford Miss. 

5. Fennimore Wise. 

6. CMckasha. Okta 

7. Waco, Texos 
e, Safford, Ariz. 
9. Reynolds, Idaho 


Records 


0.95 


10.6 


1.33 


11.9 


11.8 


0.88 


13.6 


6.25 


0.87 
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The findir^s of this analysis are sunmarized in the following: 

1) Substantial hydrologic information can be measured from 
low resolution, single-band, black and v>hite imagery 

from LANDSKT. The parameters identified and measured were; 
surface v;ater todies, land Use Type 2 and 3 surface cover 
classes, channel length, and wat'nrshed area. 

2) Of particular value are seasonal LRNDSM observations. 

For exanple, the October imagery showed the watershed when 
vegetation density was low; this made the higher- order 
streams visible, ttor.e channels could be measured frcsn the 
image than appear^ on the U.S.G.S. topogr^hic map at the 
same scale. LANDSAT imagery clearly shows that effective 
drainage density varies cyclically througtout the yeur. 

3) It became apparent that published records do not reflect 
current surface cover conditions, since the watershed's 
land use typically cilters with time. For exaitple, the 
October 1973 IANDSAlT imagery denonstrated that surface cover 
had changed fron that described in the latest published 
data (1967). 

4) Discxepancies in published data on the extent of surface 
water were easily measurable. For exanple, some inpound- 
ments had been added to the watershed under analysis since 
the U.S.G.S. map was last updated. The effect of the im- 
poundments in increasing the surface water area of the 
watershed was clearly visible. 
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1.5) Conclusions from the First Phase Effort 

1. An in^ixoved model for the prediction of peak flow events 
was structured, specifically designed to take maximum ad- 
vantage of the data and information stream available fron 
remote sensing. 

2. The itodsl was exercised to predict the peak runoff from 
nine e:!^>eriiii!ental Agricultural Research Service water- 
sheds, selected fron among a set of 158 jns^nojEnanted and 
well-described watersheds. 

3. The predictions of the new rtradel in its simplified version 
were tested against: 

a. The predictions from three of the most employed oon- 
taiporary planning models — i.e. , the Rational Pbr- 
imila method, Cook's nethad, and the Soil Conservation 
method; and, 

b. The statistical recurrence analysis of the strean- 
gage records of the nine test v;ater sheds. 

4. The results indicate that, within the range of applicabil- 
ity of its simplified version, the model spears to be an 
iaprovement over conventional hydrologic planning models. 

5. The feasibility of extracting the model inputs and para- 
meters from ranotely sensed information was identified 
Joy visual analysis of LANDSAT imagery. 
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CHMTOR II 

AfPICffiCH TO THE PfffiSE 2 EFPOKP 

The Phase 2 effort ooncentrated in tvro areas: 

1) The iitprovansnt and extended verification of the planning 
iTodel, including routing. 

2) The analysis of LfiNDSKT imagery to determine practical proce- 
dures for the €Ktraction of quantitative hydrologic informa- 
tion usable in Planning Itodels. 

The first area included four taslts. In the first task, the application 
of the peak-rate model to the test vTatershette was extended to include a 
larger, statistically significant ssrple. Estimation errors of the 
large-recurrence flow event mth respect to runoff records were deter- 
mined; the results were oorapared with estimates derived from other 
ootmonly ei^Jloyed models. 

The remaining three tas]<s aimed at the improvanent of the model's accur- 
acy by further analysis of the following important factors: 

In the second task, the tine-profile of the critical rainfall and 
runoff statistics frcan the 158 ARS watershed saiiple ware analyzed 
to ascertain regional and seasonal characteristics of peak flows 
and to determine what is the "planning rain," i.e. , the rainfall 
vAiich defines the critical outflow. 

The third task addressed the sensitivity of overland flo;^ to changes 
of physical basin parameters. This was assessed through ocxnputer' 
runs of a strip model. 
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®ie foxirhh task addressed the inclusion of '‘oon^lex*' watesrsheds. 

A routing model baaed on remotely-sensed inputs was synthesized. 

Ihe second area consisted of three tasks: 

In the first, technicmes for extraction of hydrologic data were 
analyzed. 

In the second task, an investigation was undertaken to detejonine 
the information content of each MSS band and multi-band ccmbina- 
tions. lANDSAT Imagery fron two test ^^tersheds was inteirpreted 
and ocai 5 )ared to topographic map ground truth to ascertain Which 
hydrologic features could be identified. 

Ihe third task addressed an in-depth quantitative analysis of a 
single basin fron LAMDSAT data, A Maryland watershed was selec- 
ted for this analysis because of the availability of recent, sea- 
sonal ground truth, and also because of the watershed's r^id de- 
velopment, which adds possibility of assessing the effects of 
land use changes on hydrologic properties. 

The relationships of these tasks to each other and to the overall ob- 
jectives are depicted in Figure 7, 



FIGURE 7 FLOW CHART OF CURRENT EFFORT 
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MSS RANDS 


LANDSAT IMAGERY 









CHaPTER III 


EXPaNESED VERIPICaaJION OP PEaK-RKTE MODEL AND DEVEDOPMEWT OF FOOTING MODEL 
3.1) Ejqpanded Verification of Peak-Fate Model 

Seven of the nine watersheds selected for verification in the Phase 1 
effort and which yielded the bast results were essentially '‘simple "« 
ocai^iosed of a main charmel without contributing sub-areas. An addi- 
tional twenty-four basins satisfying these criteria were selected in 
this effort. These watersheds, however, vere chosen to include a wider 
range of physiogrephic, geographic, and hydrologic characteristics. 

They ranged in area from 40 to about 2000 hectares. 

The pararaeters require by the nodel were oonputed for each watershed 
as follows; 

1) Main channel length was measured with a map v^ieel frcm U.S.G.S. 
topographic maps. 

2) The average infiltration rate was ccmputed by determining the 
weighted average of the types of soil. The U.S. Soil Oonser^ 
vation Service classifies soils in the United States with re- 
spect to infiltration rates in four classes, designated A, B, 

C, and D, respectively. Class A denotes soils with higli in- 
filtration rates even vhen thorou^ly wetted. They consist 
chiefly of deep, well-to-excessively drained sands or gravels. 
Class B denotes soils with moderate infiltration rates when 
thoroughly watted. These are chiefly mcderately deep to 
deep, and moderately well to well-drained soils with moder- 
ately fine to moderately coarse textures. Class C includes 
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soils having low infiltration rates when thoroughly wetted. 

They consist chiefly of soils containing a layer which im- 
pedes downward movaiEait of water, or soils with moderately 
fine to fine texture. Class D represents soils with very 
low infiltration rates when thoroxighly wetted. iRiese soils 
consist chiefly of clay soils, soils with a permanent high 
water table, soils with a clay layer at or near the surface 
layer and shallow soils over nearly impervious inaterial. 

An average waterslied soil class was determined by confuting 
a weighted average of the above data. It was calculated by 
assigning values of 1, 2, 3, and 4 to classes A, B, C, and 
D, respectively. BOr example, for Coshocton #194, 86% of the 
soils are type C, vdiile 14% are type B, The average class was, 
then: 

(0.86) • (3) + (0.14). (2) = 2,86 

Therefore, mean K?il class for Coshacton is approximately C, 
so the characteristics of this class were used in the model. 

Average soil class was translated to average saturated, or 
final, iiifiltration rate by using the values for each soil 
type given by H, Holtan, et, al. , for the USDAHLr-70 water- 
shed model (2). These values are presented in Table 6. 

The character of the soil vis-a-vis layers of soil which 
constitute an impediment to flow are used to determine the 
choice of the value within the range. A low value Sor im- 
peding layer of clay; a mid value for loam; and a high value 
for sand. 


mm . immmmm rm^e for soil oj^es 


Soil Class 


Final inf iltmfio?! Rate Roe^g© (cm/hr) 


I 14-76 
0.76 -.38 
0.38-12 
0 . 12 -.00 


AVMIABLE STCmGE CAPACm BY SOIL TYPE 


Soil Type 

® 7 

Available Storage Capacity m/m 

Sand 

0.29 

Sandy Loam 

0.29 

Loam 

0.25 

Clay Loam 

0.22 

Silty Clay 

0.20 

Clay 

0.1 a 
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After average final infiltration rate was calculated, avail- 
able water storage tier unit depth was corrputed. Valuefj of 
available storage capacity were assigned on the basis of 
soil type, according to Table 7, A weighted average was 
taken to determine the available storage capacity (S^) . An 
average vegetative factor (a) was oornputed for the watershed. 
First, the distribution of cover v/as determined frcan the 
data base {for example, 11% cultivated, 58% grassland, etc..). 
Then each type of cover v;as assigned a value according to 
Table 8, A weighted average was conpur.^ according to per- 
centage of each type of vegetative cover in the v/atershed. 

The ir^t set for the Holtan equation was then ccmplets. 

* 

3) Ihs average infiltration I over time was then calculated by 
the Holtan eouation; 

t .. 14 ' 

I = GI • a" ’ (S - I) * + 

3 . ^ 

4) The approximate time of concentration v’as oojnputsd by applying 
the QTpirical equation develop^ by Kerby (3) . The time of 
cancentration is the time required for rainfall from the most 
renote point of basin to reach the outlet. It therefore de- 
fines the minimum rainfall duration for the basin to reach 
p^k outflow. The Kerby formula v;as used to give an initial 
estimate of the duration of the peak rain, subject to itera- 
tive validation as described in flection 3.3, It is of the 
form; 
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Table 6 

VEGETATIVE COVER FACTORS 
HOLTAN'S EQUATION 


COVER 


0000 CONDITION 


Fallow 0.30 

Row Crops 0.20 

Small Grains 0.30 

Hay (legumes) 0.40 

Hay (sod) 0.60 

Pasture (bunch grass) 0.40 

Temporary Pasture (sod) 0.60 

Permanent Pasture (sod) 1.0 

1.0 


a) FOR 

POOR CONDITION 

0.10 
0.10 
0.20 
0.20 
0.40 
0.20 
0 .40 
0.80 

0.80 


Woods and Forests 
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c 

3/¥ 

Hheces « t±»e o£ oonoantratlon 

c 

It » distance £ecm the itost remote point in the 
basin to the diarmel, in a direction parallel 
to the slope 

s = slope 

n a retardance ooef&cientf according to Table 9. 

5) An arorage surface friction factor, Manning's "n*', was cal- 
culated jEbr each watershed. Values shoMn in Table 10 were 
wel0)ted according to the percent of ihe basin in each^cover 
class* 

A, 

6) Gwhel's Estrone Value tecimqp^ was applied to the flow re- 
cords of the e^janded sairple to oonpute the fifty-year re- 

* 

curxence flow, to be compared ultimately with predictions of 
the same event by the model. The formula is as follows: 

QgO = li + 6s 

Where: = fiflg--year recurrence outflow rate 

U = mean of flow records 

s ~ standard deviation of flow records 

6 = constant, function of the nuttOser of y^irs of 
record 

Table 11 shows values fca: Ur 6 and s for differoit recurrence 
intervals. 

7) Average basin slope was ccnpited by a wel^ted average of the 
slopes, using the data fcon A.R.S. records* 
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Table 9 

RETARDANCE COEFFICIENT - 

KERBY'S EQUATION 


Typo of Surface Value of n 


Smooth impervious surface O.OS 

* bore packed sol! ,0- 10 ' 

Poor gross, cultivated row crops or 

moderately rou^ bare surface . 0,20 

Pasture or overage grass ,0.40 

Deciduous timberland 0. €0 

Conifer timberland, deciduous 
timberland with deep forest 

litter or dense grass 0.60 
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Table 10 


REPHODUCBILITY OF THE 
ORIGINAL PAGE IS POOR 


ECOSYSTEMS 
INTERNATIONAi- IN 


MANNING'S ROUGHNESS COEFFICIENT FOR OVERLAND 
FLOW FOR VARKXS SURFACE TYPES 






Watershed Surface 

Manning's "n" 

Smooth Asphalt 

0.013 

Concrete (Trowel Finish) 

0.013 

Rough Asphalt 

0.016 

Concrete (Unfinished) 

0.017 

Smooth Earth ( Bare) 

0.018 

Firm Grovel 

0020 

Cemented Rubble Masonry 

0.025 

Pasture (Short Grass) 

0.030 

Pasture (High Grass) 

0.035 

Cultivated Area (Row Crops) 

0.035 

Cultivated Area ( Field Crops) 

0040 

Scattered Brush, Heavy Weeds 

0045 

Light Brush and Trees (Winter) 

0050 

Light Brush and Trees (Summer) 

0.060 

Dense Brush ( Winter) 

0070 

Dense Brush ( Summer) 

0.100 

Heavy Timber 

0.100 

Idle Land 

0030 

Grass Land 

0.032 








Sample 



Recurrence Interval,) 

/ears 


size, years 
of record 

10 

20 

25 

50 

75 

100 

15 

1.703 

2.410 

2.632 

3.321 

3.721 

4.005 

20 

1.625 

2.302 

2.517 

3.179 

3.563 

3836 

25 

1.575 

2.235 

2.444 

3.088 

3.463 

3.729 

30 

1.541 

2.188 

2.393 

3.026 

3.393 

3.653 

40 

1.495 

2.126 

2.326 

2.943 

3.301 

3.55^* 

50 

1.466 

2.086 

2.283 

2.889 

3.241 

3.491 

60 

1.446 

2.059 

2253 

2352 

3t200 

3.446 

70 

1430 

2038 

2.230 

2.824 

3.169 

3.413 

75 

1.423 

2029 

2220 

2.812 

3.155 

3.400 

100 

1.401 

1998 

1 

2.187 

2.770 

3.109 

3.349 
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8} !Qie $ factor Indiaatet? earlier cosipxbsd frcsn the £br^ 
inulat 

e 1 + 

l' 2 

IC3 k^ 


, kj, k^ are ratios of the overlaid to cihaoaal 

parameters of surface frictlcm factor, flow length, 
flow depth and slope, respectively. 

Ihis factor is essentially the routing ocatponent of the 
peak rate itodsl in that it accounts for channel f3.cw in^ 
tijoe of oonc^tration calculations. A detailed deriva~ 
tion can be found in Befesrence 1. 


Where, k^, k 2 


9 ) Infiltration abstractions were included in the. foim of the K 
factor vMch reduces actual precipitation to pajecipitation 
sKoess, !Ihe K term was caLculated fcy ocaiparing the 
rain rate occurring over the tuiie of concentration with the 
average infiltration rate for the sane period, ror exairple, 
the rate (P) for the 50-yaar recurrence, - duration rain 
in the Blacksburg watershed is 0, 109 mstsrs/hour. The in- 
filtration equation, using the oonstants for this watershed, 
derived as esplained previously isi 

+ 0.38 

Fbr t =3 ,52 hours {frcsn the Kirpich formula) , the infiltra- 
o 

tion rate will fall fran 7,1 axt/hr to 3.4 aivTif, with an 


Tiii:; 


PAQli; IS POOR 





• y't'':- 


f ^ 






: 


- 3 : 


avacaige value equal to approsdjnately 4.2 cr0ir* 

The K factoiTf ther^ifoirer foir this case equals: 

1 .. £ = 1 ~ ts 0,62 

P 0,109 ni/hr 

III other wards, for this particular rain ev^t, and for the 
Blacksburg wat^ahed, appcoKimately 62% of the rainfall be- 
ocmes runoff. 

10) Ml model input srsjuireraents are suimiariaed for the water- 
shed saitple in Table 12, It is clear that a great dlver- 
sil^ of hydrologic and physiogr^oihio ooiuiitions is repre- 
sented. 

11) To enable caiparison of results with the SCS model outlined 

in Peferenoe 4, an average SCS curve number wcus oonputed. 
Curve number values i^ere taken frcsn Table 13 according to 
average soil groi^. The determining factor in chooaii^ a 
curve number is vegetative c»ver, A i^ieighted average 
yielded a final curve number. 

Using the parameters developei above, the peak rate model was run for 
each of the 31 slnple watersheds. The S.C.S. formula. Cooks model and 
the Bational formula were also run for the san© sanple. The results 
are shown in Table 14: and, for the peak rate model alone, in Figure 

8. Mean errois f or the pe^ irate model were 56% compared with 62.5%, 
99*2%, and 80.3%, for the S.C.S. , COok, and Bational models, respect- 
ive!] 


OaHTJE 12 

MOEE2. JUPOP RBCK3IHSMSIPS KIR EXEESffiED I^EEERSHED SKMPIE 

( xio"^) 6ao"^) 


BASIN 


Coshocton #5 

.048 

700 

241 

Ocslwcban #10 

.049 

869 

80 

CbshDcbon #92 

.046 

2367 

209 

CJoshDcton #94 

' 

.048 

3940 

232 

ODshoctan #95 

.050 

5297 

227 

Ooshoctxm #97 

.051 

8291 

. . 

222 

Cbshccton #194 

.048 

- 

515 

163 

Ooshocfccsn #196 

.051 

902 

119 

CoshDccon #994 

.023 

17020 

242 

iCT/ra City #40 

- 

.042 

6360 

149 

Fennitnore #41 

.033 

1580 

165 

FennlnDre #42 

,035 

837 

162 

Bastings #56 

.035 

2415 

78.8 

Hastd.ngs #57 

.037 

1610 

89.3 

Hastings #58 

.0353 

1503 

75.4 

Bastings #59 

.034 

14220 

92.1 

Still\4ater #81 

.032 

1690 

107 

Wacn #83 

.035 

2400 

218 



.162 


.166 


.159 


.169 


.172 


.172 


.172 


.172 


.103 


.08 


,05 


.059 


.061 


.057 


.053 


.073 


.020 


1.410 


1.315 


1.09 


.987 


.915 


.959 


1.189 


1,321 


1,053 


1.420 


1.524 


1.535 


1,751 


1.359 


i.m 


1.916 


2.383 


1.08 


1.08 


1.08 


1.08 


1.08 


1.08 


1.08 


1.08 


1.08 


1.178 


1.139 


1.139 


1.214 


1.214 


1.214 


1.214 


1.404 



.15 

00 

■ 

.772 

.15 

.83 

.677 

.15 

.83 

.561 

.15 

.83 

.507 



.15 

.83 

.470 

.15 

.83 

.494 ; 

.15 

.83 

,713 

.15 

.83 

.612 

,15 

.83 

.680 

.151 

.79 

.495 : 

.14 

.78 

,721 

.14 

,78 

.774 

.172 

.836 

.703 

.172 

.836 

.736 

.172 

.836 

.571 

.172 

.836 

.467 

.17 

.792 

.702 

.17 

.78 

.817 


1.546 


2.646 


3.275 


5.016 


7.077 


11.96 


1.4 


i.97 


72.0 


19.4 


3.42 


2.002 


12.56 


4.91 


96.21 


6.26 


3.287 
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Table 13 

SCS CURVE NUMBERS 


Hydrologic 

Condition 


Treatment 
or practice 


Hydrologic Soil group 


Straight row 
Straight row 
Contoured 


Row crops 


Contoured 


Small grain 


Con toursd 


egumes or 
rototion 


Contoured 


Native pasture 


Woods 


61 

73 

55 

69 

49 

69 



73 

79 

55 

70 

77 








































1 





t 
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TABLE 14 






suMmsar op i®del isesults 






SEOOBDS, Q-n PEAK RATE ^SDDEa:i 

- 2 5 ERPOR 

S.C.S. 

axK 

RATIONAL 



IS ERH3R 


% ERROR 

1 



(Qeft) 

(Qcf.) 





CbahDC±on ^5 
OosliDcfeon ^^10 
Cbsbocton #92 


CoshDC±on #94 


Coshocton #95 


Ooshocton #97 


Coshocton #194 


-o^octon #196 


Coshocton #994 
lov?a City 


Pennimore %Vhl ' 


Ferminore 


I *isudLings #53-3 

1 


; Sitings #W-5 

t 


j listings #53-8 


St’. llwater 


300 tf^ 


^feoo #D 


7*62 


12,3 


5.12 


6.05 


4.63 


4.82 


10.6 


19.3 

2,69 


5.93 

1jU8 


12.98 


16.4 


10.5 


4,54 


3.57 

20,6 


13,5 


13.6 


62 

(19.91 


7 

(S.48) 


-53 

2,02 


-72 

(1,29) 


-87 
( .62) 


141 

(25.5 


-41 

(11.3) 


-96 

(. 11 ) 


-89 

(. 68 ) 

10 

(12.5) 


63 

( 21 . 2 ) 


“86 

(2.24) 


“45 

(5,75) 


-81 

(.84) 



-55 

(5,55) 


35 

(6.8S) 


-10 

(5.44) 


35 

(6.271 


-22 

(3.74) 


-58 

(4,4) 


13 

(3.04) 



76 

(22.9) 


“53 

(7.71) 


21 

(12.7) 


43 

(6.50) 


-1 

3.55) 

-47 

(10.9) 

13 

(15.2) 


68 

( 22 . 8 ) 


39 

( 10 . 6 %) 


27 

(15.6) 


112 

(10.83) 


82 

11.04) 




74 

(8,37) 


18 

( 12 . 6 ) 


-41 

(11.39) 


212 

(8.39) 


139 

(14.2) 



( 22 . 6 ) 


3 

( 10 . 2 ) 


99 

9.05) 


142 

(8.63) 



80 

(9.20) 


27 

(7.71) 


31 

( 6 . 0 ?) 


-11 

(4.29) 


66 

(17.6) 


-41 

(11.40' 


-16 

(2.26) 


-27 

(4.32) 


10 

(13.0) 


20 

(15.6) 


-40 

(9.87) 


5 

( 11 . 0 ) 


■■ 


-16 

(2.99) 


-28 

(14.8) 


3 

13.9) 


13 

(15.4) 
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TaBLB 14 (oont’d) 
StiMI®K5f OF PODBL RESlffiTS 


WISESSBSSi 

I Warn (JG 

t 

i Waco #W-2 
Waco #Y 

WacX) #Y-2 

I 

h-- 

I MbtKjuaapque #W-1 

I- . - 

Mbuqcerque 

Safford 

I 

! Safford #|i^4 

i 

L- 

I Safford #W-V 

L~ „ M-,... ■ 

Damrille #Pf~l 
Blactebiicg #1^1 
QBcSord ^10 


EBOOKDS, Q-rt PJSaK IME KDDEL 


3 ? 

m /sec./kti 


% E8K)R 
(Q50) 

-94 

(.408) 

11 

{37.J05} 

17 

(27.06) 

53 

(44.68) 


6.72 

i 55 

(10.4) 

7.61 

50 

(11.42) 

6.25 

-15 

(5.3) 

3.65 

-52 

(1.77) 

5,56 

93 

I (10.73) 

.95 

-4 

(.91) 

1.33 

1 -24 


( 1 . 01 ) 

"”9 

( 10 , 8 ) 


S.C.S. 

% ERBOR 

6B 

(7.35) 

-26 

(24.5) 

49 

(34.4) 

00 

(29.1) 

59 

(10.7) 

-2 

(7.46) 

143 

(15.2) 

95 

(7ao) 

62 

(9.02) 

125 

(2.14) 

464 

(7.3) 

-74 

(3.1) 


COOK 

RI^IOML 

% ERK)R 

i EBBOR 

(QgO^ 

®50> 

-1 

“25 

(7.35) 

(5.54) 

-26 

-19 

(24.7) 

(27.1) 

-42 

-10 


(13.3) (20.8) 


-33 

(19.5) 

8 

(7.23) 

79 

(13.6) 

-20 

(5.00) 

141 

(8.79) 

109 

( 11 . 6 ) 

478 

(S.49) 

735 

( 11 . 1 ) 

"^^29 

(8.4) 


-19 

(23.6) 

21 

(8.16) 

93 

(15.1) 

130 

(14.4) 

142 

(8.84) 

119 

( 12 . 2 ) 

405 

(4.8) 

855 

(12.7) 

-39 

(7.3) 


X = 54.1 
a = 36.5 


X = 62.5 X = 99.2 jf = 80.3 

0 = 81.8 0- =149.3 o = 163.2 


eeproducibility of the 
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‘XWo statistical tests run to detsoonins the significance of these 
results for the 31 watershed sample , The first tested whetlier the ^ ' 
variance of the paalc rate model was lower than the variances of tdie 
other three; the second tested the sigrdficance of the differences 
atong the means, Ihs results vire surmiarized in Table 15. 

In the first test, significant differences were established between 
the variance of the psaJt rate rtradeJ. and tl'e other models. 

In the secxand case, scatistical differences betoean the means of the 
psaJt rata itodel could ba shown for the Cook method, but not for the 
other tvro. The confutations are presented in i^pendix A, 

These results indicate that the reduction in variability of fore- 
cast achieved by the pealt rate model is statistically significant, 

The T'ean error of estimate vras improved also, but within the limits 
of the saitple size significesiico C6in ba demonstrated for one of the 
thr©% cases, 

3.2) Investigations to Inprove the .Accuracy of the Peak-Rate r^odel 

Davelopn^t of tlie peak irate mcdel made appairent several potential 
sources of forecast errors. It imp.row; predicirion accuracy, several 
questions ware addressed; 

i) V8hat is tils "planninc; rain,” i.e. what rainfall input comas 
ponds to the peal< flow output and what are its temporal and 


areal distributions? 


TABLE 15 


STRTISTICSai AKftLYSIS OF H313EL FOKEKSST ERKCfRS 


1) Oschran’s Itest for Bg^uaiity of Varianoa 

Test ^pothssis (B^) : Katlonal 

Accapt if G < ,7652 tg^ 


'7 

S*' test statistic (G) 


EOO Pealt Jfete 

1332.25 

- 


SCS 

6691,24 

0.83 

Inject Ho 

OKik 

22290. 4S 

0.94 

Reject Ifci 

tetional 

26634.24 

0.93 

Reject Ifo 


p% H rejected for all case's.; EOO paaJc rate variance is signi fl- 
eetly less at S5S confidence level 


2) ”T"- Test for Egnalii^ of Tteans 


Test iwothesis Uy : Meo 3 " >‘scs, Ctok, Rational 
JAcespt if T < 1.311 (tp 



X 


df 

T 



BOO 

54.1 

36.5 





SCS 

62,5 

81.8 

41 

0.52 

Accept 


Ctook 

99,2 

149.3 

34 

1, 63 

Reject 

«o 

Raticeial 

80.3 

163.2 

33 

1.09 

Accept 


" H 

BO O 

rejected for 

Coolc ircde-} 

at 30% 

confidence level 
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2) Wimt seaennal and geographic factors are jn^rtinent? Do 
rainfall and runoff exhibit propansities to occur during 
particular seasons in particular locations? What season- 
ally variant cscaiditions stiould be included in the model? 

3) What is the quantitative sensitivity of J::asin runoff to 
variations in surface pararnaters? Hoiv accurately does 
ana need to measure slope, friction factors, etc. to 
obviate the introduction of unacceptable errors. 

4) Bm should the c^ropriate basin sub-ar^ ha select^? 

Since values for each micro elemant of the basin are 
costly to measure in pra:itics, how should the i^drologic 
par-amatera be oarrbined into an average value yielding 
correct results? 

Ihs following four sections describe each question and present the 
corresponding analyses. The results from these analyses will be 
inooiporated in the model in a subsecfuant phase of this effort, 

3.2.1) Analysis of Eainfall diaracteristics 

The recurrens^ formula given in Ciiapter 1 

, a, 

X “ 

{t + d)“3 

is implicitly based v^pon toe assumption that the reiinfall is of con- 
stant intensity throughout its period of occurrence. The reason is 
that most raingages measure rainfall mass rather than rate; it is 
thus practically convenient to divide the reading of the gage by the 
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ducaticai. Equation (1) is therefore in essence an integral equa- 
tions its correct statement ^ould be the follov-ringj 



whsres 

V = rain mass, or voi.vmie, ‘.alien vn.tliin the time in- 
tei-val; t. 

In actuality, the time behavior of i uffeobs the ruiK)ff. It is 
therefore not surprising that, .in. general, no single type of rain 
event causes the peal? basin outflow. A fifty-year recurrence run- 
off, for eiiarrple, can be produced by one intense rainfall or by a 
series of lesser events. This fact, alreacfy indicated by various 
researchers, vves verified for the watershed test sample by 
ining hte relation between large- recurrence rainfalls and runoffs, 

Tha recurrence inter-ral of the runoff was calculated using tte 
Vfeibull plotting position technique, givang a regression line 
fitted through the annual runoff peaks over the period of iceoord. 
The peaks were ranJced frcm J.argest to smallest and the probabil- 
ity of the occurrence was assigned according to: 

^ n t 1 

Whsres P *= occurrence probability 

in = plotting ran}c (largest outflow - 1} 
n == years of record 


Thsss points v®re plotted as shewn in Figure 9. If tbs distribu- 
tion is assumed to be normal, the regression line can ba drawn 
assigning the mean discharge to P = 0.5. Ihs mean plus the stan- 
dard deviation, assuming normality, will have a P = 0.165, and the 
ntean flow minus the standard deviation wixl have P = 0.835, Fxxm 
this line, the recurrence interval of any discharge rate can iisi 
estimated, 

Ihe recurrenca interval of the causative rainfall was oonputed from 
tte ettpirical relation used in the rainfall module, Ekjuation (1), 
The values of the ooeffierLents in the equation were csilculated from 
records for each location, 'fhe profiles of the rainfall rate v;ere 
identified from analysis of detailed records of rainfall and runoff 
rates and irises at tin^ intervals of a minutes. In several 
watersheds, tl^ largest runoff events were identified, the enrres- 
poTKiing hydeographs charted, and the i^tographs of the generathig 
rainfall or rainfalls overlayed thereon. Figures 10a, b, depict 
typical examples. 

Oai^arisons iabween the profiles of related high-recurrence rain- 
falls and runoffs are also given in Figure 10. The oonclusion is 
that no easily-discernible relation appears to exist between the 
runoff and the mass of rainfall of equal recrurrence. Thus the 
selection of the "planning rain" is not as simple as might be ex- 
pected, Since obviously soma assumptions must ba made, it is 
valuable to Icnow how sensitive basjji discharge is to the drioice of 
rainfall profile for equal masses of rainfall. 
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As a fdjrst step, a oonputer irodel developed by D.E. Overton was 
used to gage the effect of the tsntporal profile of the rain upon 
the runoff. The model simulates the runoff resultant from a 
given rainfall ir^ut for a unit width "strip” of any length, A 
detailed derivation of the model is contained in Peferenoe 5, 

Two features are noted here: 

1) The model accounts for several important surface variables 
of watersheds. It accepts rainfall inputs of vary'ing mag- 
nitudes and tenporal distributions over strips vAxich can 
have variable slopes, dimensions, and svurface covers. 

2) The model was previously validated with accurate results 
from simulation of several hundred rainfall-runoff situa- 
tions by the U,S, Army Corps of Engineers. 

A series of oatputer simulations was itade to assess the sensitivity 
of the runoff peak to the tenporal profile of the rainfall ar«d to 
the time of occurrence of the rainfall maximun, using parameters 
typical of those of the A.R.S, v, watershed sample. Several rainfalls 
of equal volume but with the different tarporal profiles shown in 
Figure 11 were simulated. The resultant runoff hydrographs are 
shown in Figure 12. A significant differeiios in outflow: is apparent 
between the constant rain and rains of triangular shape. With re- 
ference to Figure 10, the large recurrence rain profiles taken from 
the records tend to be triangularly shaped. It thus aj:pears that 
the triangular profile should be used. 
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As regards the sensitivity to timing of the rainfall peak, Figure 12 
shows that, among the triangular rainfalls, ti^ discharge peaJc varies 
as the time of peak rain, ribiiig to a maxirnim and then falling off. 

highest discharge peak results frcan rainfalls C and D. Though 
special conditions, such as irregular basin shape, can alter the con- 
clusion, it appears valid to assume that the planning rain should 
have its maximum occur near the midpoint of its duration. 

The appropriate duration can ta estimated, as well. Figure 12 shows 
that the highest runoff rate vMdi can be expected frcm a rainfall 
of the volume input will occur at about minute 31, This gives an 
aipirical estimate of the strip watershed's time of concentration. 
Peferring to Figure 13, it can be seen that the sensitivity of dis- 
charge rate to occurrence of the peak is least in the neighborhood 
of It may be concluded, then, that the rainfall must be at 

least of duration t to ensure that the watershed reaches its ulti- 

c 

mate outflow rate. 

For the strip watershed simulated, the Kerby formula for the tirr® 
of concentration gave an estimate of of 33 minutes, satisfactory 
for setting the planning rain duration. 

Two characteristics of the planning rain are thus defined. Its tstr 
poral profile -should be triangular and its duration should be equal 
to the time of concentration of the basin. It ironains to determine 
its recurrence interval. 

Figure 14 shows hov/ this question can be addressed. It presents the 
rainfall rates plotted against their respective recurrence intervals 
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FIGURE 13 


SEr-SITIVIT-/ OF DISCHARGE jRATE TO TI^E OF DISCHASGS PEAJi 
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for the nine watershed san5>le. A marked similaril^ in the form of 
all the curves is apparent. For recurrence periods greater than 
approximately 50 years, the increase in rain rate is snail. This 
suggests that, for recurrences greater than this value, the choice 
of rainfalls is fairly insensitive. Hie planning rain for develop- 
ment of the routing model v/as therefore selected as one of trian- 
gular sh^^e, with duration equal to the time of concentration, and 
of recurrence interval of 50 years. 

3.2.2) Analysis of Seconal Factors Affecting Peak Flow Events 

Watersheds are not static, but vary with tine due to man-made alter- 
ations, and cyclically within each year because of the effect of 
seascmal changes. The response of watersheds to rainfall inputs, 
therefore, will also vary with time and season. Thus hydrologic 
planning models will be most accurate if they mirror the condi- 
tions extant in the basin during the season when the peaks are most 
likely to occur. It is therefore important to ascertain vdiether 
runoff peaks exhibit seasonal regularities and to identify the cri- 
tical seasons for the watershed under study, 

A geographically diverse sanple of watersheds was selected for this 
analysis, oonprising 15 AES basins with records longer than 15 years. 
The annual peak discharge for each basin by year and month of occurr- 
ence was recorded, and a chart prepared showing the probability that 
an armual peak would h^pen in any given month. Figure 15 exertpli- 
fies the findings, Ihe analysis permits the tentative conclusion 
that in the regions Vihere peak events are dominated by surface para- 
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inetersr (Danville, Bladusburg) the distribution is btmodal, i,e, 
flood peaks tend to occur in two distinct periods of the year, 
typically late spring and late sunnier. Watersheds in transition 
regions (Coshocton, Waco) exhibit a less marked seasonal tendency. 
The probability of occurrence of peak flows is itare equally spread 
over a six to nine month interval. Those basi is which are heavily 
sub-surface dominant (Safford, Albuquerque) shoiir single-mode dis- 
tributions. Nearly ail their peak flows occur within a three month 
period in late sunrner. 

For each saitple basin the annual paaJc discharges v/ere ranked from 
largest to smallest over the period of record. These data were 
giraphed to show larger-recurrence flows at the top, and the more 
frequent events at the bottom. Figures 16, 17 and 18 show the mo- 
dal effects described above, other conclusions which can he drawn 
are as follows: 

1) In the sub-surface dcaninant areas, flows occur during one 

season (July-Septeniber) regardless of recurrence inter- 
vals, Tlie ssirple contained five such basins with a total 
of 126 annual peaks: 86% of these occurrei in the inter- 

val above, with most of the r^iainder in June and October. 

2) Figure 19 dernonstrates an effect observed in several of 
the other watersheds. Though lower-recurrence flows 
occur in many different months, the higtest 3 or 4 flow 
events are confined to about three consecutive months. 
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Since the analysis abcxve oonsidsrad diacharga events of all recurr- 
ences, a sxibs^fuant analysis \tas undertaJcen to chart the fcehavior 


of the large recurrence flows only. Table 16 supplies the months 


of occurrence of the four largest discharges for the selected water- 
shed sainple. 


TftELE 16 FLOW HftTE RECUKI?SNCE mm. 


Safford, Arizona iW-1 
Safford, Arizona 
Safford, Arizona #W-5 
Albuquerque, M.M. , W~1 
Albuquerque, N,M. , W-3 
Ratings, N^raska W-3 
Hastings, K^nraska W-5 
OcshDcton, O, , §5 
*Ooshocton, O. , iJlO 
*waoo, Itexas, C 
*Waoo, Texas, W-2 
*Waco, Texas, Y-2 
•Stillvreitex, CSda., W-4 
ICMa Cii^, Iowa 
Pennlitore, VJisoansin, W-4 


Those waterst^ls which stowed no 


1 

2 

3 

4 

Sept, 

Aug. 

Aug. 

July 

Aug. 

Sept, 

Aug, 

Sept. 

July 

Aug. 

Aug, 

July 

Aug. 

Sept, 

Aug. 

June 

Aug. 

Sept, 

Aug. 

July 

July 

May 

July 


July 

June 

June 

June 

June 

Aug. 

June 

July 

June 

Sept. 

July 

April 

Iferch 

Aug. 

April 

June 

May 

June 

J^iril 

Itov. 

ifey 

June 

April 

S^jt. 

April 

Oct. 

July 

May 

July 

July 

June 

Aug, 

Aug. 

Aug. 

Aug, 

July 


seasonal trend are starred. The 


Coshocton, Waco, and Stillwater basins, as noted ^rlier, are in 
transition regions and shov; no seasonal sketmess, Fran tto data 
above it would appear that, accept for transitional regions, peak 
flovre tend to occur during a j;articula£ season. Therefore, the 
watershed conditions extant in that season should be reflected in 
the hydrolcgic model. 


This hypothesis was tested further by performing a similar analysis 
for discharge volume rather tlian rate. Table 17 gives the xEjnth of 


occsunrence of the maxiinura two tour discharge volvme. J\gain, the 
font hipest events for the period of record are presented. 


TPSm n Blig'g K>H]MB REiCUgRBCB BRNK 



1 

2 

3 

4 

Saffiard, Ariaona, W-1 

Sept, 

Aug, 

Aug. 

Sept. 

Safford, Ariaona, W-4 

Aug, 


Aug. 

Sept. 

Safford, Ariaona, 

July 

Aug, 

Ai^. 

Aug, 

Albugcerque, M.M, , W-1 

Aug, 

Sept. 

Aug. 

Aug. 

Albuquerque, N.M, , Vh3 

Sept. 

. Aug. 

June 

Aug. 

Hastings, Nebraska, Vh3 

fte/ 

June 

July 

July 

Hastings, KebrasJ<a,W-5 

June 

June 

July 


Qjahocban, 0. 

June 

Aug. 

June 

March 

*C3ostoctan, 0, ^)10 

June 

Sept, 

July 

i^ril 

Texas, C 


Aug, 

J^il 

Dec. 

Ofexas, W-2 

Nay 

/^ril 

June 

Kov. 

*Waco, Texas, Y-*2 

May 

April 

March 

June 

*Stiliwater, C2d.a,, W-4 


MSy. 

July 

Oct. 

Iowa City, leva 

Jid.y 

Aug. 

July 

July 

*FenniinDre^ Wisoonsin, W-4 

Aug. 

July 

July 

Jan. 


Ths Starred watersheds once more e^diibit a weator seasonal tenj3sncy. 
Witli the exception of ths Fentiiirore basm, ths rate and voltEte sam- 
ples are tlie sane, Oortparison of the rate and volume table yields 
ths following intormaticms 

1) Of the 60 sairple points, 58% of ths rate and volums teaks 
occur in the sane month, and 80% occur within one moniii 
of each other, 

2) Only seven of the samples did not occur in the same season. 


©lereibre, except in transitional regions, flood rates and volumes 
@{hibit a propensity to occur in a particular thr^-nonth season, 
typically spring or late sunrosr, A planning model, then, will be 
mcona liJcely to produce acceptable estimates if the condition pre- 
sent in these seasons are reflected. 
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3,2.3) Analysis of Sensitivity of Runoff to Surface Parazi^ers 

purpose of this emlysis is is-rofold? 

1) to dstamins to/ accurately surface characteristics must fce 
measured for input to models, i.e. , what are aco^table 
errors in estimation of values of hydrologic paramstsrs. 

2) to develop a rationale for ''averaging*' surface parameters - 
for sumring tlie values for each point on the watershed into 
earn conputationally manageable unit. 

In the first phase of this effort, sensitivity to sub-surface drivers 
was analyzed via analog csxrpjter Efficgrams. This ta^c, therefore, c^- 
ter^ K>ley upon investigation of sensitivity to stirfac^ parameters. 
The strip model described earlier was employed, because it provided 
the capability to vary all surface drivers and detennina ch^ges in 
the resid-tant ruraoff . Computer runs ware rasde to assess the sensi- 
tivity of discharge to slopa and surface friction. 

All cxsrputec runs assume constant rain. The ocEnbinatic^n of triangu- 
lar rain profile and varying surface Esrameters is reserved for a 
subsequent phase of this effort. Also, the oonputer runs assumed a 
strip of fixed total lergth and a f ijted duration and intensity of 
rain, Thus the results presented hereinafter are to be oonsider^ 
as indicative only. The go:isa:ali 2 ation of the results to strips of 
aridtrary length, and rains of arbitrary intensity and duration is 
reserved for later phases of this effort. 
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1, ghe Effect of Slopa on Runoff Peak; 30 runs 

h series of ccjinputer nois vsara made to deteimina the effect 
of various slopes -topical of rural watersheds (0.01 - 0.20) 
on both a single-plane strip and a t».t»-plane strip. Tiie 
first series assume! a strip of unit vrldth, constant sur- 
face friction (n = .05) and lengtii (1 = 30Sm) and variable 
slope, A constant radnfall (P - 8ar\/hr) was mathonaticaily 
applied to the strip and the runoff hydrogiraph recorded fer 
slope values of .01^ ,05, .15 and ,20, The results s5:s3W 
that the value of slope iicpacts the timing of the occurr- 
ence of the discharge peaJc, i,e . , the time of ocncentra- 
tion. As expected, hi^er slopes produce a storter t^; 
this, in turn affects the design duration of the planning 
• rainfall. Figure 20 shows, for exat^sle, that the 205 slope 
produces a time of <x>ncentration of ^proximately 16 ndn- 
xrtes where a 5% slope gives a t of 20 minutes. At a given 
recum^os intenral, say 50 years, and fcr a given regional 
rainfall fonnulation, say Oxford, Mississippi, this change 
in slope results in an approximate decrease in peals: rate of 
11%. For the sst® set of assimptions, however, a diange in 
slope from 5% the seraa proportion to 1,25% girnas a 38% 
decrease in p>eak rate, Tia cxaiclusion is that slop>e must 
ba m^sursd more critically for flatter basins (S < 5%) . 

In addition, fciie results (Figure 20) sliov? that the duration 
of the p^k is a function of slope, asynptotically approa- 
ching 20 minutes. As shown in Figure 21, for slopes less 


f. 



L8 


24 


30 






















gressively Increasing nutiber of streaght planes. Figures 


28 and 29 present the results. 


1)16 approximation of the watershed profile by a single plane 


produced errors of in the peak rate duration. Itie inclir- 


sion of a second • _ s greatly in c recused accuracies. Further 


approximation did not, however, produce significant improve- 


ment, This suggests that only major changes in watershed to- 


pography need to be aooounted for. anall fluctuations (<5%) 


in slope can be smoothed over without loss of accuracy. 


2. The Effect of Surface FTictio > on Runoff Peak; 


The surface runoff models uses Mamirg’s "n" as a meMure 
of flow resistaj^oe due to surface roughness. Overland flew 
was modelled for both the one and two plane cases with 
varying s>irfaoe rou^uiess. In the single plane run, the 
strip chosen was identical to that used in the slope sen- 
sitivity runs. All paremeters were held constant (slope ® 
.05; 1 = 305m; P = Sonvhr) while Manning's "n" was varied, 
(n * ,01 to .10) r^aresenting sxirfaoe covers from ooncrete 
to heavy forest. Surface friction was found to affect 
peak discharge rate as shewn in Figure 30, The time of 
oanesntration chosen for the planning rain will vary frra 
9 CGraph A) to 39 minutes (Graph E) dependi n g on the 
Mannlrjg’s "n" selected. Errors in discharge resultant 
fr o m inaccurate surface friction estimation will be most 
severe for low values of Planning's "n". It may be oon- 
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cluded then, that the nore developed a watershed is, the 
more critical is the quantification of surface roughness. 
Further, peak duration were found to depend on surface 
roughness as shown in Figure 31, Peak duration is sensi- 
tive to surface i huess principally for more resistive 
surfaces {n > 0.04). For rural watersheds, a realistic 
range of Manning's "n" is from 0.03 to 0.10. In this re- 
gion, a 25% estimation error will result on the average in 
about a 50% inaccuracy in runoff duration. 

In the twor-plane runs, the single strip was again divided 
at the mic^int. Both sectioas were given identical slopes 
and leigths, and various conbinations of "n's" were 
modelled. The steacfy uniform rain was applied. Figures 32 
through 37 show the sensitivity of the two plane model to 
surface rou^ess. Insults similar to those of the single 
plane sirttilation were found. The insert of surface friction 
measuranent upon estimation of t^ is most critical for low 
values of Manning's "n". 

To date, oonputer analysis has been conducted for a basin 
strip of average dinension and for a typical rainfall input. 
The principal findings were; 

a) Duration of runoff peak is most sensitive to slope at 
values less than about 10%. i*7hen slope of watersheds 
are greater than 10%, an average value can be approxi- 
mated, In flat basins, though, more detailed ground 
truth should be consulted. The same sensitivity 
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applies to estimates of time of concentration *“ tie 
critical region consists of tlie lower slopes, 

b) Irregular watershed slopes can be s^sproximatad by 
straight planes where slope variations are less 
than about 5S; otherwise, each strip should be di- 
vided with the outflow of tiis i:pp3rnost bscatving 
the input to the lower, 

c) rurxjff paalc rate is irost sensitive to surface 

friction in basins vri-th lower resistivities (n < 
0.04) , vMle the duration of the peak is most sen- 
sitive to M^er friction. It follows that a re- 
motoly-sensed estimate of surface cover needs to fce 
adequate enough to separate it into classes with 
similar values of jyiatnning's "n". Referring again 
to Table 11, it bsccmas patent that a rarote sen- 
sor should be able to separate forests from fields, 
fields froin soil, and soil from urban areas, for 
ejsnple. "\ 


._^From the above, the following criteria can be established for the 
selection of itodel parameters ; 

1} The rainfall ir^t should be of tnang^ar sliapa, have dijora- 

\ 

tion equal to the tim-a of concentration, and ba of recurr- '• 
encs equal to fifty years. 


i 
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2) The model parameters should reflect the coiiditions extant 
in the basin during the flood-probable season. 

3) The slope module should be particularly accurate for flat 
basins; the surface friction factor shou.ld be differen- 
tiable among classes with large differsices in Manning’s 
"n". 

3,3) DevelopBent of Routing Mcdule 

In large watersheds, or in those crrnposei of several tributary streams, 
the assunption that all areas of the basin contribute to the outflow 
l^fdix>graph simultansously leads to errors. The hydrograph of each sub- 
t-^atershed can differ from those of its nei^ibors in teamporal distribu- 
tion, magnitude and dixration. The overall outfit frcsa the basin is 
the ocartoinatian of these hydrographs, appropriately arld^ tr> account 
for the time lag required for runoff from eu-ch to reach the basin out- 
let, Figure 38 illustrates the differences betv>sen temple and ooinplex 
basins. 

C3omples watersheds rypically contain than c./, prediominant channel. 

The outflcvt> frrjm eacti sub-are i does not drain directly into the main 
channel; rather, some flows first to seoc'ndary streams. The contribu- 
tion of certair sub-areas is delayed. Figure 39 ,’llustrates tiie 
effect. An effective h/drologic model accomplishes mathsnatically 
what has taen done graphically in the Figure: it ocHi^rates the hirdro- 

graph of each sub-watershed or sub-aren and sums them according to a 
tine-delay (routing) function. The routi.ng^ model can provide the 



FIGURE 38 


ILLUSTRATION OF SIMPLE AND COMPLEX WATERSHEDS 


COMPLEX WATERSHED 
NEAR MUSCODA, WISC 


SIMPLE WATERSHED 
NEAR MUSCODA, WISC 




■ 













Chi 


( Hf bn>n< 


m 


m 


mm 


FIGURE 39 


SUMMING OF lELAffiD HYDRDGRAPHS 


TOTAL DISCHARGE 



94- 


water resource planner a tire picttore of the buildup of the peak 
flow event, vdiereas the sirrple model generally yields only its 
maniinum rate. 

Tte approach to developing the routing nodule of the peak-rate itEx^ 
was to aipiOKiiiate a watershed by a series of unit strips and then 
to sOTi the hydrogr^hs according to a time delay function. ^Ehis tech- 
nique was graphically summarized in Figure 39. Ihe model, therefore, 
consists of an overland flow caipon^t (unit strips) and a channel 
flew ccrpcaient (lag time function) . In line with the intent of this 
effort, the rodel was designed to keep the oenputing hardware to a 
xninimuai, 


Ohe model appcosdmates overland flow on a unit-width strip and outputs 
the discharge hydrograph. The overland flow oorponent assunnes outflcw 
to be a- ftmetion of (tepth at the discharge point: 

Q= a hj' (1) 

3 

T''here: Q = outflow rate, 1 /time 

= depth of outfall, 1 

m,a - constants 

It also holds that the outflow at steac^ stats will equal the rain 
rate times the area: 

Q = iL * 1 (2) 

tflhere; • i = effective rainfall rate, accounting for 

ilfiltration losses 

(width = 1) 

oh^^ = il>; or 


Therefcare: 


(3) 
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h = 

o a 


( 4 ) 


0316 average depth over the strip will equal the sunraations of the 
f tepths at each pDint divided by the total length: 


h = 


J 


o a 


(5) 


adhere: 


h » average depth 
£ t= length betajeen 0 and L 


* (-) 
'a 




1/bh 


at 


( 6 ) 


* If the release of inmoff at the outfall adl^res to the Maaming rer- 
lationf then iti = 5/3; 


h= (i) 


vs 




i 1 » 

= 's> 


8/5 


(5/8) 


(7) 

( 8 ) 


c. n 3/5 • 

I s 5/8 h 

8 a ' c 


(9) 


h _ 5 


( 10 ) 


The storage of water on the strip equals the average height times the 


area: 


S a hL = 


( 11 ) 


Where: 


S s= storage 
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lUffetrentiatdiig botii sides yields: 


dS _ 
St ~ 


dh - 
O #S rx 

sr t 


( 12 ) 


The change in storage an the strip vdll also equal the infla^r rate 
minus the outficrrf: 

Inflcw e (rain rate) x (area) = iL (13) 


Outflcw » (tfetffling relate) » 


1.49 


5/3 


(14) 


n 


Iflhere: S - slqpa > , 

n - feiction factor 

InfLoi? - Outflow (IS) 

UU 

Etxxa equation (12} ; 



ds 
St ^ 

ST 

( ^ Xi) = 

s Inflow - 

Outflow 

(16) 

Therefore; 

51. 


=s iL — 

1.49 

it 

^5/3 

(17) 

or, 





^5/3 




*^o 

- 8 4 

- 2.38 

/s 


(18) 


dt ’ 



n 

IH^ 




Ecpiations (14) and (18) define the overland flow ccsnponent of the model. 
It was progranmed to rm in an iterative fashion on the lazger-c^city 
land calohiators (Hewlett-Packard HP-65) . A sanple oui5>ut ^pears'’ in 
Pigar© 40.. "5Sie required inputs are slope, surface friction factor, 
flew length, and- rainfall rate. 23a first three are potentially measur- 
^e frem remote sen^^rs. 
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A 3X)Uting function was then added so that the xunoff oonttihfation of 
sevHcal strips can be accumolated and a watershed modelled in totality* 


^Ehis was acccnplished ty introducing a lag time for- each strip depenr 
ding t^pon its ^stance frcan the basin <?utlet, the slope and friction 
of idle channels, and tte rain rate. The particular equation etsdoyed 
was tak^ from Overton ^ferencs 5) and is of the form: 


tag Time (Tj^) 



0:63 



(19) 


Where; i = rain rate, ft./sec, 

n = Manning's "n" for the channel 
L « length feon outlet to strip, ft. 
s = channel slope, ft,/fh, 

The outpit of a strip with length fcom the outlet of 100 feet, sur- 


face friction factor of .05, slope Of .05 and rain ra^ of .00007 feet 
per second will have its lag time ccnputed as follows: 

5 / (.05) (.00007) (1000) 

8 (.00007) ^ 1.49 /T^ 



r 


s 506.1 secs. 


= 8.4 minutes 


®ie outflcw hjidrograpih from this strip, as ocitpufced fcy the overland 
flow ccBpcnent, will be delved by 8.4 minutes and then added to the 
discharge of the remainder of the basin to yield the overall basin hy^ 
drograph, Figure 41 graphically depicts this modelling process. The 
ccnplGte form of the routing model tterefiore; 

•. 1.49 

n - O 


n 


( 14 ) 



FIGURE 41 Lagging Strip Hydographs 
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*^0 _ 8 . 2»38 /T“ h 

F ^ — 

T « 5 ± nit' 

^ 8 1.49 /s 


( 18 ) 

(19) 


output of equation (18) will Ise tiie depth of flow at the outlet at 
time » t. value is inp\ 2 t to equation (14) to ccnpute outflow 

rater nlso at titpe- - t. 03ie l^drogr^ih w^l be formed £com the Q's 
for all t’s. Xts'< contrihution to total basin discharge will he deter- 
minedr in turur by equation (19) . 

Previously r sifa-surfaoe sei^itivity had been e}^mined with the 'analog 
‘coippater and surface sensitivil^ through the siitple strip nodel. A 
final of sensitivii^ required a complete rouinng model for its an- 
alysis. The sensitivity of discharge to ptysiogxcphic parameters/ 
namely basin sh^pe and areal distribution of rainfeU/ was detenmined 
thoooui^ the use of the newly developed model. 


A simulation was run to quanti^ the effect of watershed shape tpon 
peak discharge rates. Five hypothetical basins of 400 hectare area 
were modelled/ each having a different diape. Figure 42 shows the 
shapes and their cojODespanding discharge ouiputs for a constant rainr 
fall input. As expected, those basins with large percentages of their 
total aarea oLose to the ovtlet poxaduced higgler discharges. The sh^Je 
factor in this ejcanple accounted for ip to a^ 20% difference in outflow 


peak. 


The time to peak rate was also affected, though not significantly. 
These findings confirm oiar earlier assertijon that watershed shape is 
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a driver of peak flow events. Since shape is a discharge driver, a 
planning nodel should be capable of accounting for it. The routing 
model provides the Cc()ability. The watershed will be simulated by a 
number of strips viiich can have different lengths (as well as slope 
and Manning's "n") , so any shape can be approximated. 

A second simulation was performed to ascertain the sensitivity of 
outflow to the areal extent of rainfall. To this end, a square basin 
of 400 hectares vras assumed, ard different rainfedJ. rate ooiribinations 
applied over its area. A hi«^ (ISoVhr.) and low (7,5an/hr. ) rate 
were simulated for a 30 minutes duration. The basin was divided in 
thirds and the high rain was effectively moved from the outlet to the 
ridge line, Figure 43 gives the results. 

Mgher peaks accrued frcm harder rains nearer the outlet. This is the 
same effect observed in the shape simulation in that higher peaks re- 
sult v;hen a greater percentage of the rain can reach the outlet quickly. 
Further, it is possible to get a dual rather than single peak resultant 
feon the position of the rainfall relative to the outlet as shown in 
the Figure. Uiough the particular oatibination of rain rates and posi- 
tions vdrLch produce the maximum outflow will vary with basin size and 
shape, vSien basins to be measured are very large or when rainfall char- 
acteristically occurs only over small areas the effects of the "wDtet 
case" stoiild be accounted for. 

Again the routing model is particid.arly applicable to modelling vari- 
able rain rate. Each unit strip can be given a separate rainfall in- 
put. Where detailed rainfall network data exists, then, it can be 
modelled very closely, 




^Chef £oUcK«lng sunnarizes the findings and results of the develops 
ment of the planning model: 

1) Ihe peak rate model has been further verified by ext en s io n 
of the test sairple to 31 watersheds, The results of the 
model found to be a significant improvanent over si- 
milar (xsRnonly-used models. Most of the model parameters 
can potentially be measured from ranote sensing stations - 
slope, surface cover, and physiographic measuranents. 

2) The analysis of rainfall characteristics and their seasonal 
modifi^irs has allowed conclusions. First, thou^ no 
direct relation exists between the "x” year rain and the 
discharge of like recurrence, oiough evidence exists to per- 
mit estimation of a *'plannir^ rain,” The rationale for its 
s^ection has been presented. These rains formed the inputs 
to the peak rate model and produced reasonable results. 

Second, it will behoove the water aresouroe modeller to con- 
sider the seasonal perturbations to rainfall-runoff phenom- 
ena, Individual geograEiiic regions will have varying peak 
flow characteristics. These factors provide valuable clues 
to oonstruction of the planning model data set. Satellite 
ranote sensing offers an inprovenent over ccaiventional 
sources of ground tsnith in this regard, Shcmld it be deter- 
mined that the model should reflect late suBmer conditions 
in a basin, for ^lanple, lANDSAT imagery from that period 
can be acquired and analyzed. The most readily available 

Or' TO. 
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source of aacial pihotograptos, cm the other banfl^ "sees" a 
pgarticular area only once every several years* 

A routing \co6b1 has been synthesized for jiodelling of “oont' 
pleac” basins. It offers principal features. Firstf a 
great deal of flexibili ty is provided for approximation of 
vjatersheds V7ith internally variable lydrologic parameters. 
Ihe nuitber and pl^sical characteristics of the strips can 
be tailored to closely approxiitiate conditions extant in the 
■watershed to be modelled. Second, anple opportunity for 
raraotely-sensed input is provided. Slope, surface friction, 
strip and channel dimensions are potentially measurable 
frcsn ranote sensing stations, . The model has beai success- 
fully applied -to assesanait of the sensitivity of snanoff 
to basin sh^ie and areal distribution of rainfall. 
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BSDEQiaSIC 

This ta^ included three principal criipanents: 

1) Analysis of tihe state of the art of tednigues for measuring 
surface characteristics of hyftrologic significance and of 
their costr inplsientation time^ egudprosnt and skill require- 
ments. 

2} Evaluation of the specific hydrologic information content 
of the four lAKDSAT hands and determination of \iMch bands 
or cczihiinatdons of ban^ are best suited to measuririg each 
modal ii^ut. 

3} Quantitative hydrologic analysis of a ocnplete watershed 
using lANDSKT imagery and available ground truth. 

4.1) Analysis of tine State of the Art of tleasuraiient of Hydrologic 
Paran^ters frcm Paootely-Sensed Il&ta 

Thie investigation focused rpon measuraoent ir^chods for those reciotely 
sensible hydrologio factors of most imtediate value to watershed model- 
ling* 

1) PhysiogrEphic Basin Maasuranent 

2) Surface Oover Identification & Classification 

3) Soils Classification 

The anal^is could carried to significant depth thariks to the simul- 
taneous availability of an effort sponsored by Goddard Space Flight 
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Oentei:; Dr. Vincent Salcnonsonr l!^5''20567, "User Bequirements and 
User jUsceptance of Current and Next*K3eneration Satellite Mission 
and Sensor Ccnplenent# Oriented Ttaward the Mjinitoring of Watet Re- 
sources," of DecanbeTf 1975 r with which the burden of this research 
was shared, 

4,1,1) Sfeasurement of Physiographic Basin Pasrameters frcau Esnotely 
Sensed Data 

Various physiographic nBasuarements which serve as direct inputs to hy- 
drologic models are possible ftxan lassiDS®!!? imagery: the measuranant 

of watershed area, overland flow lengths, drainage density; and the 
estimation of channel dimensions* 

Watershed areas are d^ined fcy the ric^e line, i,e, , the contour withr 
in ^diich all rainfall drains to the same outflow point. In regions 
vhere relief is pronounced, iddge lines can be directly discerned frcm 
aerial photocp^phy or satellite imagery. Where slopes aare irst so steep, 
topographic and slope m^s are valuable interpretive aids. With the 
boundary delineated, area can be measured with a planimeter or calibra- 
ted grid. Plow lengths can also be directly measured. The watershed's 
sh^pe, vhich has been shown in this effort to be one of the drivers of 
discharge, can also be discerned. 

Pigure 44 illustrates the method for the definition of the tflatershed 
area. The t^jpermost set of images danonstaates how the ridge lines 
can be identified frtsn pronounced relief. Shown is a western California 
watershed of 72,320 hectares. The left-hand image shows that a signifi- 
cant contrast exists between the noarthsastem (light) and southwestern 
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(dark) slopes o£ the ste^ ridges ocsnprising the basin boundary. In 
the right-hand image, a contour is drawn around the ridge line der 
fined by the reflectance difference. The area can then be measured 
by overlaying a calibrated grid and counting the divisions vnLthin the 
line. 

The lovser pair of images shows an a3:ea vhere slopes are less than 10%, 
so that rriief is not apparent. What is apparent, however, are sev- 
eral streams, each with associated drainage structure. The streams 
appear dark in the Eand 5 image while the land shows tp light* By 
following the streams visually it can be determined which ones drain 
to the same point, thereby appaxadmately defining the watershed. The 
boundasy location can be refined by interpolating the distance be- 
tween streams which drain away from each other. The line on the ric^t 
hand lower image was drawn by dividing the distance between streams in 
half. This appcoxdmation of the ridge line resulted in area measuare- 
ment errors of order 5% or less. 

It is possible to determine piysiographic data fam laiJDSftT imagery 
to an accuracy ccmparable to that available frctn mediurarscale maps, 
iigure 45 shews a section Oi an U.S.G.S. lt250,000 scale topographic 
leap of an test watershed at (3iidkasha, CKkla. The elevation con- 
touars give the altitude of points arxjund the central stream. Ety • 
connecting the points of highest elevation, the watershed boundary 
can be obtained. The area of the basin was measured by planimeter 
fran the map to be 54,857 hectares. 




Comparison of LANDSAT Imagery and Topographic Map for 
Determination of Watershed Boundary and Area 
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lo«<@r illustraticoi is a image enlarged four times to match 

the scale of the topograEihia map* 2his vjatershed is relatively flatr 
so ridge lines could not be relied i^on to define area. The boundary 
was instead determined by using the drainage divide and vegetation 
Shade changes as guides. She watershed area was mensurated by a polar- 
planimefcer to be S6r470f a 3% overestimate with respect to the topogra- 
phic map. 

'i 

Should the apprcpriate dues not Ise present/ the satellite image can 
be projected onto a topographic map and aligned using identifiade 
features CcoadS/ 3akes) as landmarks. Ohe ridge line deriv)^ from 
the topo m^ can Ihen be transferred to the image. 

The length of channels affects the ratio of overland to channel flow. 
Their particsular arrangement within the watershed ^temdnes the 
amount of time that precipitation spends running over the land versus 
in the channel/ and therefore the lag time involved in ta^insmitting 
flow to the basin outlet. The length of the Channels is a variable 
entering the caiputation of basin time of concentraticn in certain 
l^^drologic modds. In others/ drainage pat-hem aijd density are ccm- 
pared with basin area to develop time delay histogram for overland 
flow/ because it determines the lag time Sar the contribution of each 
watershed sub-area to reach the outflow point. Watersheds with 
higher drainage densities have more Channel length per unit area and/ 
therefore, water will travel a shorter distance overland before rea- 
ching a strean* Since channel flew is generally faster than surface 
flew, a waterslted with hi^ drainage density will reach its outflow 
peak (jiickar than one of equal area with a loror ratio. 
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Figum 46 detiDnsbrates the dit^rbance of drainage density in h^dro* 
logic phenomena. Ihe upper image is a 400 ^uare mile section of 
Band 5 LBNDSSS? image, photographically enlarged ^proximately 2.5 
times. It shews a number of watersheds with a high ratio of stream 
length to rmit srea^ Streams appear light and land surface daadc. 

A r^resentative drainage density for the area is appmsdmately 1 
Mlcmeter per square Mlometer, Hhe SKact value fbr each hasin can 
be measured by first determinii^ basin' area as earlier described 
and then caiputing stream l^igth within the area with a map wheel 
or similar measuring device. 

The lower image is an area of like dimensions £ccm a higher and 
dri^ location. It is apparent that mqch more surJ^ce area per unit 
stream length eidats in this watershed than in the previous image. 
Stream length was measured by map wheel frem the image to be 210 
kilemeters* Drainage density here is approximately 1 kilcmeter per 
5 square kilometers. 

Figure 47 deononstrates the utility of the several IAND6AT 3?ands for 
measuring channel parameters. The images show a Minnesota watershed 
as seen in each band, To measure total stream ler^[th, best results 
are achieved using Bands 4 and 5, because higher ord^ (smaller) 
streams are visible vMch cannot be seen in the other bands. On the 
other hand, to measure channel width, Bands 6 and 7 are preferable 
because the contrast between the stream and the surrounding vegeta- 
tion is much inpreved. The main channel is particularly visible, 
spearing light amongst darker vegetation. 
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4«1*2) ffeasurenmt and Classification of Surface Cover fron Barotely 
Sensed 3^gecy 

.« 

Tlie and dlstribirfcion of siarface cover influeaace the friction en- 
countered h |7 precipitation eKcess viiile running overland* Soil t^pSf 
another surface paraoneter, has laore iirpact t^!on suihsurface processes. 
Vegetaticrx and land lase can be classified remotely in sufficient de- 
tail to permit the assigrsnent of quantiiative hydrologic values. 

®his will be demonstrated later in the'3?eporb. The classification of 
soils is more difficult, however, because in m^y cases the soils 
•cannot be seen due to vegetation cova:. 

S^drologic planning models typically divide watersheds into zones vdth 
similar vegetative characteristics, or, klternativ^y average the vege- 
tation factors across a catcJment and utilize a ’’lumped" parameter. 
MorOTver, v«?getation class and density are used as indicators of water 
retention and infiltration. 

Figure 48 indicates the ability of imagery to determine diff- 

eient vegetative cover on the watershed as ccnpared with other ground 
truth sources. The ipper photo indicates a section of 1$24,Q00 (ori- 
ginal scalel US(^ typographic quadrangle including the fiRS test waters 
shed in Black^urg, Virginia. In the lower aerial photo (original 
scale 1:45,000} the field structure and hydrologic land use is readily 
apparent and with limited grovmd truth can be successfully interpreted. 

Field stnmture can be likewise interpreted from the I^SM3SST image. 
Ifote particularly tl^ odd-shaped field near the center of the aerial 
and satellite photos. Classification fron LSilDSAT data will generally 
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require seme ground trutdi^ such as aerial imagecy, Sox correlating 
tbs surface cover class and r^lectance level. Figure 49 illustrar 
tes the procedure jBor incorposrating the I^MIsaaj-derived data into 
ccaiventional hydrologic inodels. Once the water^ed is partitioned 
into hydrologically distinct classes such as ra^ crops, fallow, 
forest, impervious etc* , the hydrologic parameters can 

he ascribed to the individual subwatersheds. 

^e technique is illustrated for a hypothstical watershed for three 
important paramsbers; Manning’s "n”, iMch is a measure of the imr 
pedance of the watershed to overland flow; Holtan's *'a'*, which is a 

t ^ 

measure of the ilfiltacation potential of the t^tershsd; and the SCS 
curve number, vhich is directly used to determine the fraction of a 
given rainfall event which appears as direct rtnoff * 

4»1)3} Bemote Classification of Soils 

Hydrologic models igpically group soil into homogeneous hydrologic 
classes of different permeabilities, porosities, etc, Areal extent 
of soil type is used to determine infiltration rate and potential 
moisture content. Soil pemeabilii^, porosity, and oonductivil^ 
serve as inputs to ihe infiltration and ev^xjtranspiration oenpo*- 
nents of the models. 

Soils classification must in many cases rely upon inference rather 
than direct measurem^t since tte soils arre often obscured by vege- 
tative cover, . Soil association must typically be inferred from 
knowledge of the surficial vegetation, Sesne interpretation is pos- 
sible directly frem the Imagery since soil wetness affects photo- 


ASSSGNMENT OF HYDROLOQSC PARAWIETERS TO 
VEGETATION CLASS 



VALUES OF HYDROLOGIC PARAMETERS 


MANNING'S "N" 

HOLTAN'S "A" 

SCS CURVE NUMBER 

SMOOTH ASPHALT 

0.013 

FALLOW 

0.10-0.30 

ROW CROP 

CONCRETE (TROWEL FINISH) 

0.013 

RAW CROPS 

0.10-0.20 

(STRAIGHT) 67-91 

ROUGH ASPHALT 

0.016 

SM. GRAINS 

0.20-0.30i 

(CONTOUR) 65-88 

CONCRETE (UNFINISHED) 

0.017 

HAY (LEGUMES) 

0.20-0.40 

SM. GRAIN 

SMOOTH EARTH (BARE) 

0.018 

HAY (SOD) 

0.40-0.60 

(STRAIGHT) 65-88 

rlRM GRAVEL 

0.020 

PASTURE 


(CONTOUR) 61-84 

CEMENTED RUBBLE MASONRY 

0.025 

ibunch grass) 


LEGUMES OR 

PASTURE 

0.030-0.035 

TEMPORARY 


ROTATION 

CULTIVATED AREA 

0.03S-0.040 

PASTURE (SOD) 

0.40-0.60 

NATIVE PAS- , 

SCATTERED BRUSH, HEAVY WEEDS 

0.045 

PERMANENT 

0.80-1.0 

TORE OR 39-84 

LIGHT BRUSH & TREES 

0.050-0.060 

PASTURE (SOD) 


RANGE 

DENSE BRUSH 

0.070-0.100 

WOODS & FORESTS 

0.80-1.0 

WOODS 25-79 

HEAVY TIMBER 

0.100 



• 

IDLE LAND 

0.030 




GRASSLAND 

0.032 
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graphic darkness in MSS Band 7. Ancillary data such as vegetation 
type, geogr^ohic location, slope and pisndmity to water bodies fa- 
cilitates classification and msnsuration» 

Pigtire 50 sha^s a section of the SCS Soil Surv^ for Marshall County, 
Mississippi. The surv^ provides a map of general ssil associations 
and an acre-fcy-acre couniy^-wide soil classification. Tte figure de- 
picts the soil association map: it alone is generally marginal for 

hydrologic mooelling since hydrologic parameters vary widely within 
a soil association. 

The lower figure is a I.ANDSKP Band 7 image of the same area takai 
August 1973. The differences in shade allow broad soil association 
classification. A clear oorrespondencje exists betv^en ireflectance 
in the image and soil association as shown in the SCS maps. Addi- 
tional detail useful in finer classification, is however not readily 
apparent. Thecefore, in models rat specifically designed for remote 
sensing iipput, additional ground truth, such as the soil survey, is 
required to assign values to subsurface ii^ts. 

4,2) Specific ideologic Analyses of LANTSAT Iiragejy 

4,2,1) General Principles 

The ofc^ective of this task was to determine the eirtent to which in- 
formation directly applicable to hydrologic models can be gleaned 
from the satellite data. The task oonsisted of an assessment of the 
information content of the MSS bands and of the quantitative l^'dro- 
logic analysis of selected test basins, ^secifically; 
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FIGURE 50 


Photo interpretation of Remote Sensing’ Imagery 
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for Soil Classification 


Soil Conservation Service Soil 
Association Map 
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1) tibe identificatdcai of hydrologic features of watersheds, 
nanjely, stream courses, forest cover, surface water and 
iiri^jermeable areas. 

2) The asaigment of two quantitative parameters used in the 
routing model for the two basins - Mannings* "n" and 
Holtan*s *'a'*. 

^Ehe effort focused upon esctraction of the pertinent information by 
relatively single, low-cost, visual methods, enploying black and 
\4iite Landsat imagery, 2hou^ more sophisticated techniques such 
as ooDoputer-aided analysis have yielded good results, lOOsystems* 
faaaHiarity with anal! users indicates that the benefits from re- 
motely^sensed data vjill be realized earlier if sinple yet adequate 
in&rmaticm-extraction procedures can be identified, 

^Ehe previous ta^ indica-tad that at least a minimal level of ground 
truth greatly facilitates the interpretation of ranote imagery, be- 
cause it assists in establistog a correlation between surface ob- 
servables ar»i reflectance levels. 5hie degree of availdaility of 
such ground truth data is rot the same everysdiere, OOpographic m^>s, 
soil surveys, and a^ial photography can be cibtain^ for large areas 
of the U.S. Often, however, all three do not exist simultaneausly 
for particular watersheds. In spite of this mevenness of coverage, 
even limited ground truth is valuable. Therefore, the capability 
of siperinposing Landsat imagery onto ground tanath data at the same 
scale is of major lirportance to successful extra<3tion of hydrologic 
parameters* 
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Ob acconplish this function, an image viewing system was assatibled, 
capable of projecting single or multi-band IMDSPIT imagery onto a 
work surface. The satellite imagery was magnified up to forty times 
and could be traced or overlayed ground truth maps or photo- 
gr^is. The system developed is a low-oost, highly reliable visual 
analysis tool of the type readily accessible to practicing hydrolo- 
gists. 

4.2.2) Analysis of Hydrologic Information Content of I'JiNDSKP Bands 

Two of the JfflS ej^ierimental watersheds were selected for visual pho- 
tointerpretation of LANDSAT images. Those chosen - the Thorne Creek 
basin near Blacksburg, Virginia, and Watershed W-10 near Oxford, 
Mississippi - were singled out due to the availability of adequate 
ground truth in the ferra of topogrcphic maps, soil surveys and aer- 
ial photography. These aids were applied to verify the interpreta- 
tion from LKNDSAT imagery. 

The same methodology was used for both basins. Using the projec- 
tion device, the Imagery in each band was overlayed on a tppograjiuc 
map which the basin boundary had been drawn. The LANDSAT im- 
age was matched in scale, position and orientation to the map using 
proninfint physical features (roads, rivers, etc.) as referaace 
points. Once matched, selected portions of the LANDSAT image were 
visually ocnpared to the ground truth to establish "^training sam- 
ples" of reflectance correspondences. Next, the ground truth was 
removed, and the photointarpreters attoupted to identify all sur- 
face features and to measure those requiring quantification. 
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4.2*2, 1) toalysis of the Bladcsbitcg Watesrshed 

Figure 51 reproduces an aerial photograph of the HLaciksburg water- 
shed at the scale of approxiinately Is 60,000, The basin is 1235 
hectares in area; it contains primarily agricultural fields, plus 
small stands of forest. Urbanisation is insignificant. The fol“ 
lowing describes the Jesuits of the photointerpretation in each 

IfiNDSM? band, and compares the detail discernible in lANDSAT imagery 
to that available in other sources. 

Band 7 (0.8 “ 1,1 yen) 

a. Surface Water “ The Blacksburg watershed csontains no standing 

water areas of significant size. The LSNDSflT scene in vMch 
it appears does contain a large ('vl,000* width) river: this 

was readily discernible and used to align the imagery with 
tite ground truth map. 

b. gbreamcourses >- Without the topo map, identification of small 
streams was exbronely difficult. When the IiMDSaff image was 
overlayiad on the map, some stream patterns became more ^par- 
ent “ approximately 20% of the streams in the \-)atershed were 
detectable. 

Vegetation - Vegeta'fcive detail is not .readily detectable in 
Band 7. Agricultural fields could not be delineated, ^bst 
vegetation appeared a uniform gray shade. Correspondence 
with topo map information was fair. 55 hectares of land 
were identified as forest, whereas 71 are reported in the 
topo map. 




Impeoneable Areas - The only in^jermeahle areas readily iderr 
tifiable were roads. Primary four lane highways appeared 
(^kr secondary roads were partially visible. Sane railroad 
right of way could be seen, i^roximately 25% of total high 
way miles could be seen. 


Band 6 C0.7 ~ 0.8 um) 


a. Surface Water - The large river in the scene was visible, 
though not quite as dark as in Band 7. This river again was 
used to align the imagery. 


Streamoourses - No streams could be discerned within the 
watershed without the aid of the topographic map. When 
overlayed, however, approximately 20% of total strean length 
could be idaitified, appearing a slightly darker gray than 
their surroundings. 


Ve getation - 5he more heavily, vegetated areas of the water- 
shed appeared darker in Band 6, though forest could be iden- 
tified positively only with the aid of the topo map. Forest 
area was estimated by planimeter at 56 hectares frcrn the iitt- 
ages and 71 from the topo itBp. 


d. Boadoourse - With the topo map overlayed, approximately 35% 
of total highway length could be measured. R»ds appeared 
dark in Band 6? ti»y could also be identified by shape. 




Band 5 (0»6 ~ 0,7 ym) 

a. Surface Water - The large river was visible in Band 5: csai- 
trast was not as good as in Bands 6 or 7, Visibility was 
sufficient to use the riv^ for alignment, however, 

Strearocourse - With the image projected on the map, more 
streamcxAirse became visible in Band 5 than fran the p3?ece~ 
ding two bands. 7^>prcodiiately 22% could be measured vis- 
ually. 

Vegetation - Band 5 appears to offer the most vegetative de- 
tail, Color ranges frm dark to very light and sane field 
' shapes are distinguishable, I%>rest area was measured by 
planimeter as 52 hectares, 

Raadoourse - Hoad detail was minimal without the ground 
truth. Contrast is not as good as it is in areas where 
roads are cut through forests. When overlayed, though, the 
image revealed sectior^ of pciinary road, approximately 
5490 meters of the 12,505 meters present could be discerned. 

Band 4 (0.5 - 0.6 tim) 

a. Surface Water - Surface water is obscured in Band 4. The 
large river was only faintly visible and could not be used 
to locate the watershed. 

b. Streamoourse - No streams could be measured even witli the 
tqpo map bb refer^ice. 


c 


• Veggfaation •* Vegetative detail is significantly worse than 
in Band 5. Msst cover classes appear an alinost unifom 
gray shade^ with a few fields having hi^ reflectance. One 
such field could be located and correlated to a ground truth 
aerial photo, fbrest acreage was not measurable, 

d. Boadoourse - Ctnly occasional sections of road could be seen 
and only in ccatMnation with the topo map. Thoi:^ two sec- 
tions were used to align the map/ no significant quantita- 
tive measurecients could be talcen, 

4. 2. 2. 2} Analysis of the Oxfcrd Watershed 

Figure 52 shows the Oxford, Mississippi test basin. Ihe watershed con- 
sists of 2227 hectares in north central Mississippi. The land is low- 
lying, characterised by a preponderance of forest cover with seme open 
fields and no rohanisation. 

Band 7 {.0,8 - 1.1 urn) 

a. Surface Water - According to the 1970 USGS maps and the 
ASCS aerial photography, the warershsd contains 105 ponds. 

Of these, 22 are larger than i acre in area. In lAMDSAT 
Band 7, 11 ponds were visible, ranging frem 1.25 to 8.0 acres. 
49% of the total standing water area was visible. The visi- 
bility of the water bodies was not directly related to the 
size of the water bodies, i.e. , seme ponds of area 1 to 2 
acres were clearly discernible while other larger ones were 
nob. Distinctively ^aped pcaids were successfully used to 
align the lANDSAT imagery and the topo m^. 

S.S5^R0DtTC^1T..p^ 
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Strearssourse Streams per se were not visible# but the 
streamcourse could be identified £rcm the differences in 
vegetative cover associated with the water route* Wide 
areas around streams appeared light, but did not follow 
the tcjtal laigth of the streams* It was possible to mea- 
sure 10,370 alters of streamoourse out of a total of 

I 

20,283 meters, drainage densil^ w^ not measurable on 
Band 7, Streamoourse contrast was valuable in aligning 
the Band 7 imagery to the nep. 

Vegetation - Coarse vegetation detail was visible cri Band 
7. It was determined from aerial, pbotograpihs that tte 
forest cover appeared a medium gray vMle less heavily 
vegetated areas were much lighter. Only tx*) gray shades 
. could be easily separated by visual analysis* The oacres- 
pondence batfween forested areas on tie IffiKDSKP images and 
tee map was generally good but not exact. T!te borders be- 
tween forests and non-farests Hatched in^iarfectly. Forests 
crctprised 1730 hectares on the LHHEJSKT images and 1321 hec- 
tares oh the map. 

d. " Impermeable Jhreas - Ihe only inpermeable areas clearly visi- 
ble were tee largo: roadbeds. Visibility of roadways in 
Band 7 was minimal - no quantitative neasures could be made. 

Band 6 C0*7 - Q»8 tm) 

a, Surface Water - Some ponds wsre visible on Band 6 hut con- 
trast was poorer. Six ponds were identified ranging frciti 
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1.25 to 8.0 aares* These six vjeare a subset of thuse i<3entl- 
fled Band 1 , i.e. / five of the inteznediate size ponds 
viere obsqured. These six ponds represented 41% of total sur- 
face water area, 

b, Stoearocourae - ^i^iearance of stieamcourses on Band 6 was 
siinilar to Band 1 , but the 'amount of detail visible was re- 
duced, General ccjntrast was slightly poorer than that of 
Band 7, 8450 matecs out of 20283 maters total strearooourse 

vere measurable. The streams thaniselves were not visible, 
though a canal approcmately 50 f^t wide was identified. 

In Band 5 the streamcourses and ponds were u^d to align 
the image. 

Vegetation - Forest and non-forested areas \^ere visible, 
although agreement vdth the map was iniierfect. The heavily 
vegetated areas were basically a single shade of gray. Ti^ 
detail was slightly tatter than in Band 7. Forest area was 
estimated from Band 6 as 1704 hectares, a 29% discrepancy 
with respect to the USGS ne^. 

d. Impermeable arejs - Koadcourses were far more apparent in 
Band 6 than in Band 7, The rocuis in the scene are of secoir 
dary and teitiary size. The roads thennselves v?ere not visi- 
ble but the ocmtrast bstvi^en forest and clearing around the 
road was apparent, i^proximately 50% of the ridge lin^ 
corresponding to the path of the roads was visible, 7320 
met^s of a total of 12810 meters of roads were identified* 
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Band 5 (0.6 0,7 m) 

Surface Water *» Surface water was invisible on Band 6# ifo 
ponds could be identified. 

5treamcou3:3e - alie streamcourse in Band 5 was presented less 
confca^st than in either Band 6 or Band 7, It was nore diffi- 
cult to locate the basin since streams and ponds were net 
satisfactory guides. ®3ough the oontxast was severly reduced, 
more strearacoiarse was visible than in Band 6 and Band 7, once 
the Imagery had been registered with the map. I^d 5 poovided 
an inproved estimate of total channel length over either bands 
6 or 7. 

c. Vegetation - Subtle differences in vegetative cover were 
^jparent, Ebrests were darker than fields, but greater dif- 
ferentiation was possible within the forest class. In this 
respect more detail is present than is given in the map, and 
the correspondence of forest area is greater. 

d, Impermeable iUreas - Foadoouraes have hi^ reflectance in 
Band 5 and contrast is marked. Tertiary roads are visible to 
a great esetent, and were used to align the watershed imagery. 
Without the roads, location of the watershed would have been 
esfcranely difficult. J^jproadmately 80% of the ridge-line 
corresponding to roadoourses was disoemible. 

Band 4 C0.5 - 0.6 um) 

a. Suarrace WSter - Nb standing water was visible in Band 4. 
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b. Sfareamoaiirse •• Sti'earncourses were lights than fbrest but 
contrast v/as not as sharp as in Banots 6 and 7. About 6% 
of total str^nicourss x^as detectable in Band 4. 

Vggetation - Overall contrast is poorer than Bands 5, 6, 
or 7, light and dark correapondenoe with groxand truth was 
similar to Bands 6 & 7, but visibility xtfas severly inhi- 
bited, Ibrest area, therefore, was not nieasured in Band 4. 

d, Itnpgm^ble Areas - l&me but not all secondary and tertiary 
roads were visible in Band 4. Ibtal length visible ti®s 
visible vias 10,d80 met^s, Boadcourses appear light and 
were used to align the imagery, but contrast is much worse 
than in Band 5. About 30% of the ridge line ccaareaponding 
to roadcourse was dear. 

In additional analysis, the photointerpreter took a single surface 
characteristic - forested area - and attempted to quantify it in 
bands 5 and 7 and in two multi-band caiibinatioins, Tito forested areas 
xflere napped and subsequently oonpared to ground truth data {USGS topo- 
gr^hic m^) to chedc dassification accuracies. Figures 53a through 

6 d^ict maps of forested areas constructed from Landsat imagery by 
utilising the reflectance levels of the bands. It is dear that Band 

7 greatly overestimates forest area. This is consistent vsith the 
single band analyses p>resented earlier. Band 5 can differentiate 
more detail anong surface covers and therefore yields a better esti- 
mate. Satisfactory results (15,1% invoitory error) also came frcm 
the intersection of tonds 5 and 7. 
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FIGURE 53 
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Table 18 presents the emission and cxmiission errors calculated 
jEcom the forest area estimates. In the invaitory node (errors cal- 
culated by taking the difference between observed and actual areas 
and dividing the result by actual area) , Band 5 yielded an error of 
only 1%. 

4,2, 2.3) Oonclvtsions from the Analyses 

Table 19 suntnarises the findings made on the Qoeford and Blacksburg 
watersheds. The following oandusions can be made regarding the re- 
latlwi value of each bai^ limitedr of course, to the &;ampl^ ana- 
l^ed. 

1) .Band 7 appears to be best for identification and measursment 
of surface water area. This is because of the very loir; rer 
flectance of staaading water in the 0,8 *• 1,1 imi range and 
ocaisequent high contrast with its surroundings. Though no 
significant urbanization odsts in either basin, analysis of 
the renairder of the USHCSiaT scene showed Band 7 to be good 
also for identification of urtan land use. 

2) Band 5 offers much more information about vegetation than the 
other bands. Whace vegetative cover typically e^pears only 
as one or two shades of gr^ in other bands, Band 5 often 
yields twice that nunher. 

3) The information derivable frott Band 6 is correlated with that 
of Band 7; likewise, Band 4 is correlated with Band 5. In 
both cases, the detail in the former has proven inferior. 

Band 4, hovrever, was found useful in the measurement of road- 


oouroes. 


TABEE 18 


Itotal Area 
Measured 

Qnoisslon 


]iiventoiy 
Mode Error 

'Actual area 


FOREgr/VBSEaaTIVE AREA C(WARISON 


Band 7 


1730 Ha. 


6 . 1 % 

(81 Ha) 


31 . 1 ^ 
(411 Ha) 


315« 


« 1296 Ha. 


Band 5 


1309 Ha. 


9.3SS 
(123 Ha) 


21 . 2 ^ 

(281 Ha) 


1 . 01 % 


Intersection 
5 & 7 


1123 Ha. 


2 . 8 ^ 

(37 Ha) 


18.95! 
(250 Ha) 


15. 


1911 Ha. 


2.1SS 

(28 Ha) 


38.3?! 
(506 Ha) 


44.5?! 


( fim U.S.G.S. 1:24,000 Scale Otopographic 


IXXCB: The U.S.G.S, is not necessarily accurate. Thus the 
results qiDted should be ocaisidered as indicative, but 
not conclusive 
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j !: 'ABLS 19 I CO^g ARISON OF MEASUBED mUROLOGIC SURFACE CHARACTERISTICS TO USGS TOPO, MAPS 


STREAMCOURSE FOREST ROADCOURSE SURFACE WATER 

VISIBLE (ft,) VISIBLE (Ac.) VISIBLE (ft.) VISIBLE (Ac.) 


i 

I 

i 

1 

Blacksburg 

Oxford . 

Blacksburg 

Oxford 

Blacksburg 

Oxford 

SB 

Oxford 

J 

t 

iBAND 4 

i 

* 

|B 


1792 

* 

36,000 

m 

0 

1 

( 

ii:AND 4 

14,000 

36,000 

129 

3232 

18,000 

12,000 


0 

BAND 6 

12,000 

28,000 

138 

4210 

14,000 

24,000 

mm 

33 

BAND 7 

12,000 

34,000 

135 

4273 

10,000 

0 

0 

39 

MAP 

58,000 

163,000 

175 

3265 

41,000 

42,000 

B 

80 


TaSIE 20 COMPARISON OF HYDROLOGIC PARAMETERS ESTIMATED FROM LAHDSAT TO THOSE 

BSTIMAa’ED FROM GROUND TRUTH DATA 




Holtan* s "a" 

Maiming’ s ’’n” 


Va, 

Miss. 

Va. 

Mi-ss. 

BAND 4 

Ik 

.59 

■sk 

.053 

BAND 5 

,42 

.73 

.033 

.071 

BAND 6 

.43 

.86 

,033 

.084 

BAND 7 

.43 

CO 

.033 

.084 

BAND 5 & 7 

- 

.70 

- 

.065 

MAP 

.44 

,76 

.034 

.072 

A. R« S • 
RECORDS- 

,42 

,56 

.039 

.060 


* No quantitative measurement possible 


































final- step of the anal;^is was aimed at detenninlJig the accuracy 

of the measura^^ats of Manning’s "n” and Holtan’s "a" from LSNtSOT 

C 

uiegery* TSiese two parameters were first oca^juted frcm the M!S re- 
cords, The two coefficients were calculated fron tlie information &ar 
trached £con each band. For the Bladc^Durg watershed, the best es-ti- 
raate came fccm Band 5, though it was only marginally better than the 
vali:es obtained from the other bands, 'All estimates are satisfactory 
for mod^dlng purposes. The best single band for the Oxford watershed 
was Band 4: the Band 5 and 7 intersection. was an iiigorovanent over 
other single-band figures, Besults are givioi in Table 20. 

The overall conclusion is that since eacii band provides an information 
increment, the sum of the information present in each should be ^jplied 
to determining ttiQ model’s parameters. This can be acccirplisl:^ by an- 
alysis of eacli band individually, as w^ done for the Oxford and Blacks- 
burg basins, or by cxnpositing the information into a ’’synergistic" im- 
age, as will be described in the next section, 

4,3) Quantitative I^drologic Analysis of lASIDSST Bnagery 

The 3 r^ults of the preliminary visual analysis of the Black^mirg and 
Oxford watersheds were sufficiently premising as to warrant its ex- 
•tension to a detailed study of a -third basin, possessing hi^-quali-ty 
recent ground truth. The objective was to determine how many surface 
features could be identified and measured; to assess -the accuracies 
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adbievable in areal nieasuretnent in the inventory and land use rradesj 
and to ascertain the hydrologic information content of ocstpssite imr 
agery, The Muddy Etcanch Creek in Montgomery County, Maryland, 
selected for this study, for the following reasons: 

1) A set of multi-tenporal ground truth is available in the form 
of highraltitude color infrared aerial jdiotographs of high 
quality. These were obtained through the cooperation of the 
Iferyland Department of State Planning, The photography was 
taken within two months of the LRNDSKP pass analyzed. 

2) Soil Surv^s and 1:24,000 scale topographic maps are availr 
able, 

31 The waterslted is gaged. 

4) The watershed is sufficiently large - approjiimately 5000 hec- 
tares •" to contain a varieiy of land uses, thus presenting a 
good test of the remote sensing capability to identify a di- 
VBsrse set of surface covers, 

5) Other local studies of this watershed for other purposes are 
available: thus a good inventory exists of physiographic and 
lydrologic data. 

6) The watershed is rapidly urbanizing and subject to frequent 
flooding, resulting in hi.ch roonstary losses, 

7) The 'basin is sufficiently close to permit detailed on site 
visual examination if jrtjuired. 


Figure 54 shows a map of tl^ Branch watershed, surface cover 

contains fields and meadows, seme of which are cultivated, and sane fallow, 
plus other se^al land uses. The areas jjimediately adjacent to the 
main sti^am are predcniinantly forested. The northeastern end of the 
water^ed contains part of the city of Gathersbu3rg and therefore, re- 
sidential and industrial land uses. Several lakes are present, ranging 
fron oner-half to approximately 5 hectares in ai?ea. Figure 55 is an 
aerial photograph showing the location of these suirface covers. 

The iaro bands containing the most hydrologic information, MSS 5 and 7, 
were first analyzed individually: the accuracies of measurement of im- 
portant hydrologic parameters frem Landsat imagery were determined by 
conparison with the accuracy achievable frem aerial photography. ®ie 
anal^is technique was visual. Thirty-five millimeter sections con- 
taining the Mud^ Brandi Watershed were cut out of LMSlDSKr 9" X 9" im- 
ages, suitably mounted and projected by means of the device described 
earlier. '’Training'' oorrespciidences betweai surface cover and visible 
color were established? the basin was then classified using only the 
landsat imager^'. The findings of the photointerpreters are presented 
foUrarlng. 

Band 7 CO. 8 - 1.1 ym) Slack & White 

Vegetation detail was not good. CJontrast among vegetated areas was low, 

Forests and fields were a rtediuni shade of gray mtii foresees appearinc: 

only slightly darker. Thus, forested area was generally overestimaiied , 

# 

similar to what had ocemnred for the Oxfoird basin. The photointerpre- 
ter measured the forest area by planimeter to be 2591 hectares , as ocm- 





Cornell 





pared with 1519 measured frcm the aerial photographs. Band 7 is mudi 
better suitei for identi:^ing and measuring surface water. Forty one 
lal^ are present in the watesrshed. Eleven could be seen easily; 
two mojre were discjemible when the topographic map was overlayed. 

Uiese thirteen, caopri se 67% of the total surface water. Roac3s were 
not particularly visible on Band 7, R major interstate highway which 
traverses the watershed could be seen, as could a second major artery: 
no others were visible. Streamoourse visibility was minimal and could 
be identified only with the aid of the topographic map. When over- 
Ictyed, 5490 meters of streamoourse out of a total of 79,300 wsre dis- 
cerned. 


Band 5 C0.6 - 0.7 pm) Blade & White 

The river banks of the Potomac River vMch ^jpears within the scene 
containing Mudety Branch Cre^ vrere not sharp. Xn^iermeable areas - 
standing water and urbanization were also depressed. No lakes or 
ponds could be seen. Only the large interstate highway was visible. 
Conversely, vegetation detail was substantial. It ^'^3S possible to 
more accurately differentiate bebveen fields and forested areas. Area 
was measured by planimeter to be 1703 hectares , a oveirestimate 
with respect to ground truth. Differences in vegetative shade oould 
be v^ed to delineate certain other physical features of the watershed. 
Power line cuts ootald be easily seen; as could the location of a 
building situated in a snail cleared area. Streamcourses per se oould 
not be discerned, but were apparent insofar as they were correlated to 
forested .ureas. 
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Slngle~band analysis o£ Madc^ Branch thus produced results congruent 
with the earlier findings; Band 7 is nest useful for identification 
of surface water and urban uses while Band 5 is nost suited to vege- 
tation classification, 

Multi-band Imaqery Analysis 

Multi-band images were prepared for anal^^is using a diazo processor. 

9" X 9" black and white images were transferred to a sheet of 
transparent film and eicoded in shades of a single color* A Band 5 
image/ for example , wras reproduced in shades of magenta. Band 4 and 
7 diazo images were produced in yellow and cyan, irespectively. This 
color cenibination places the individual images 120 degrees apart on 
the *'o61or vdieel" and thereby heightens the contrast and readabiliiy 
of the oonp^site. 

The individual transparencies of tl^ three bands v©re then oombined 
into a layered sandwich. lEhe images were registered using the re- 
ference marks available on the Landsat imagery. Tests with cali- 
hrated test sheets demonstrated registration errors to be signifi- 
cantly less than a pixel. The advantages of thus producing diazo 
cen^osites were found to be: 

1) The information contained in all the spectral bands is 
liade visible simultaneously. 

2) It is well-known that the eye is more amenable to separa- 

0 

tion of color than of shades of gray. In fact, the addi- 
tion of color provides a marked inpiovanent in visual 
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analysis over black and white iitages. The three primary hues 
conibine in the ocwposite to produce a large number of distin- 
guishable colors* 

3) The diazo films can be developed and oombinBd to stjoess diff- 
erent colors and shadings. This permits the interpreter to 
*'tune” the ccnposite images to accent the contrast of desired 
features. For exanple^ a ocotpOsite with a dark magaita ccm- 
ponent in Band 5 obscares vegetation oontrast, but heightens 
surface water visibility. 

4t) The process is very inexpensive, A four-color oonoposite can 
, be generated for a total cost of approximately $1,00. 

A slide-size section containing the watershed was cut out from the com- 
posite and projected as previously described. 

The first ocotposite analyzed was made fron LSNDRAT images taken in 
October, 1973, The ground truth aerophotos were taken in Decsiber, 
1973 at a scale of 1:130,000. Five distinct land tees could be sep- 
erated; forests, fields, lakes, bare soil, and urban areas. The 
LAllDSAP oortposite was visually interpreted and thematic maps pre- 
pared for each cover class. Description of the analysis and accur- 
acies achieved follows, 

FOHBSTED AREA 

All vegetation appears red or red-orange; forests display a darker 
3red shade. In the inventory node, forest area was measured as 1242 
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hectares fern the LaUDSKP oonposite versus 1292 from the aerial photo. 
This estimate of forested area is significantly Isetter than that ob- 
tained from single-band inagery. 

FIEEDS 

On the aerial photograph, cropped and non-cropped fields were dis- 
criminable; in the Lf3SIDSK2 ocnposite tl^ agricultural areas were too 
similar to ordinary meadows to be separable from these. Fields could 
usually be s^>arated easily from forests due to their lighter red- 
oraxige ^lade. Scsie of the boundaries between fields and forests or 
bare- soil were unclear and had to be estimated. Area of fields was 
measured as 2852 hectares fran the LANDSAT imagery, a 6% underesti- 
mate ootpared witii aerial photoground t ath. 

LAKES 

Standing surface water showed the best contrast with its surroundings 
and was therefore the most easily identifiable surface element. Lakes 
down to an area of approximately one acare were visible, though all 
lakes extant in the basin could not be identified from the LSNDSAT com- 
posite. Surface water appeared almost black in the ooirposite and was 
typically surrounded by orange fields, illustrating the value of con- 
trast for identification and mensuration. Percent total surface water 
area was overestimated slightly from the LANDSAT imagery diK to the 
dcaninance of border pixels in the low Band 7 reflectance of water. 


^^myuojprrrr^ ^ 

r nrp, 


-147- 


BaRB SOIL 

In the aerial photograph, land without vegetative cover showed tp white. 
In the color ccn^josite, areas identified as hare soil presented 

two distinct colorations, ftost areas ^ipeared light blue: sctne dis- 

j£Layed a naich darker blue colcaation. The reason for tliis difference 
is not oainpletely clear; it spears that moisture content could be 
causing the effect. Both li^t and dark blue areas were summed in con^ 
puting bare soil area. 575 hectares acres were measured frcm the conr* 
posite versus 547 frcsn the aerial. 

uBBawaREas 

Two distinct shades were visible for developed urban areas. Sane lo- 
cations appeared pitMsh while others were bluish. Both types were 
less uniform in color (i.e, , appeared mettled) than any other cover 
class; this feictor allowed their separation and identification. From 
the aerial photo, the type of urbanisation (residential versus non-re- 
sidential) could be seen and a correlation was attenpted bett^en land 
use and color. Though the arelation did not always hold, there appeared 
to exist a correspondence between color intensity of development. 
Xiarge-lot residential areas contai n ^'egetation (and therefoire appear 
pinkish) than more crowded sites possessing a less amount of natural 
surface cover. Total urban area was measured to be 557 hectares frem 
the LSNDSM! composite and 649 hectares from the aerial. 

Figuce 56 is a black and white reproduction of the October LS45DSKP 
composite. The 2cwer graphic presents a summation of the thematic 
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maps jiKToduced for each cover class, Ihe areal esctent of each land use 
was measured by ovsrlayiiig a calibrated grid and counting tha squares 
incd.uded in the particular cover class. The number of squares was then 
converted to hectares and errors calculated. Soma areas had earlier 
been estimated with a planimeter, however, it was found that the grid 
technique produced comroensurate accuracies while reducing analysis 
time. The results obtained in the inventory mode are reported in Table 
21. Th^'were calculated ty subtracting actual area from measured area 
and dividing by actual area. The errors shown, therefore, represent 
differences in areal measurensnt rather than in location. Table 22 pre- 
sents the results in the land use mode. As expected, these oorrespon- 
ding errors (area correctly classified minus actual area, divided fay 
actual area) are higter than the inventory errors. However, in most hy- 
d3X)logio models, the parameter lased is the percent of a watershed in 
each surface cover class: the inventory mode errors apply in this case. 

Table 23 presents a breakdcfwn of errors by category. LaJces, for exam- 
ple, were most oft^ confused with fields. This was not due to their 
having similar reflectance characteristics, but rather to the fact that 
the border between the two was unclear. The lakes supplied very low 
reflectance and therefore dcminated the border pixels. This led to 
confusion with the surrounding fields and consequent overestimation, 
Bbrest, field, urban and soil oambinations exhibited similar b^vior. 
Though their color sha53es in the conposites were easily distinguishable, 
the boimdaries were not as sharp as in the aerial photography. For ex- 
ample, two forested areas separated by a small tract of fields tend to 
be classified as all forest. 
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LfiND USE MOEE ERRORS - MUDDY BRANCH BASIN 


gAHU E 22 

i 

LAND USE MODE 


Urban 

Area 

COMMISSION 

ACRES 

OMMISSION 

ACRES 

LAND USE 
ERROR 

156 

472 

- 24% 

Forest 

614 

778 

- 23% 

Lakes 

22 

16 

- 28% 

Bare 

I 



Soil 

684 

546 

- 46% 

Fields 

1,088 

1 

i 

1,548 

- 20% 
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Ihe p63?cent areal ccnsposition of each class of cxwer of the Muddy 
Branch watershed was oonpxited frcsn the aerial jihotogr^iiiy and from the 

LSNlsar analysis. Table 21 ^ows that the largest difference between 
the two was one peroaat. 


Manning's ”n" and Holtan's "a" were again calculated on the basis of 
L?iNDSflT and ground taruth data. Table 24 depicts the findings: 


Manning's "n" 
Holtan's "a" * 


TABIg 24 

LiMJESfiT 

O3MP0SrEE 

0.037 

0.489 

for non-inrban 


SERIM- GROUND 
TRUTH 

0.037 

0.480 

only 


The color ocst^site clearly gave results equal to those of the aerial 
ground truth. The accuracies achievable from LSNDSHT in the assign- 
ment of hyd3X>logic parameters were found to be quite satisfactory. 


The quantitative analysis of LfiNDS?^ imagery for Muddy Branch water- 
shai demcnstrates that accuracies sufficient for hydrologic modelling 
Ci be obtained through relatively sinple visual means. It was fur- 
th . determined ffcm the investigation of the Cscford and Blacksburg 
watershed that two factors can iirprove classification results: 

1) The quality of the data used as ground truth is important. 
The aerial photogrcphs used for the Muddy Branch analysis 
were moi® recent, and therefore more reliable, than was the 
map data applied in the Oxford and Blacksbiirg investigation. 
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Msresver, aerial photogr^jhy gives a nicture of vAiat is 
can the surface rather than a graphic representation, thus 
facilitates cxaiparison to satellite imagery. 

2) The application of the several MSS bands iandividually and 
in ccsiposite form pcovides far more discximinability and 
therefore better accuracies than does the use of a single 
band. 
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CHfiFEER V 

SU^3^EIQ^ OF FIKDINGS J^ID CDIKXUSIONS 

Shree taaJcs were undertaken in this effort: 1) the validation of the 
peak-rate nctdel on an expanded set of water steds, 2) the developnait 
of a routing Jiodel fear oanplex basins, and 3) the quantitative hydro- 
logic analysis of IiWJSSff imagery. IJhe findings and results are de- 
scribed in the foUcwli^. 

5.1) Expanded Validation of the Peak-rate Mndel 

A thirty-one watea^ed san^de was selected with significant geographic 
and hydrologic diversity, Ihe prediction svpplied ty the irodel was 
tested on. each, and the output ccnpared to the records and to forecasts 
ocnputed hy using three other conventional plaruaing models. The remote 
sensing model gave inproved variability and accuracies ccmnsnsurate to 
the other three models, Kfean errors for the raoote sensing model were 
approximately 50%, The renote-sensing model in its current h^lemen- 
tation applies to "sittple" basins - corposed of a single predominant 
channel, and devoid of significant sub-basirs. 

AdditicQ'ally, some potential sources of modelling error were identi- 
fied and therefore a nunher of pertinent questions were addressed. 
First, the "plannit^ rain" had to be defined. Subsequent analyses 
led to the conclusion that this rain could be best approximated by one 
of triangular shape, having a duration approximately equal to the time 

of ooncentration of the basin and having a recurrence of ^roximately 
¥ 

fifty-years. Second, seasonal characteristics of peak fla^ phenomena 
were investigated to ascertain what impacts they mi^t have on a model. 


It was disoovered that different geogr^jhic regions ejdiibit varying 
seasonal properties^ but that, within a region the characteristics 
are sijtdlar, Ihose basins located in subsurface-dcrainated areas, 
for example, shoi^ a propensity to produce peak discharge in a two to 
three month period in late suniner hence requiring increased satell~ 
ite coverage during this period. She model, therefore, should mea-? 
sure the pl^siographic (drainage density) and hydrologic (surface 
cover, soil moisture, etc.) conditions which exist in the critical 
season. Finally, sensitivities of surface parameters were examined. 
It was shown that the runoff rates were sensitive to slope primaarily 
at low slopes. Further, it was found that sensitivity to sxarface 
function requires that a remote sensor be able to classi:^ surface 
cover into categories with similar values of Manning’s ‘'n”. 

5,2) DeveloEitiant of a Routing Model 

The need for a model to ti:eat "oGnplex” basins was identified above. 
A model was developed which approximates the v;atea:shed by a series 
of strips, each having its own set of surface and rainfall para- 
metets. The output of these strip® is sunmed using a simple time 
delay function which accounts for the length of overland flew and 
the hydrologic characteristics of the channel. The ooiplete model 
was ^iplied to analysis of the sensitivity of runoff to basin slope 
and areal extent of rainfall. It was discovered that both are sig^ 
nificant and should be provided for by the planning model. The 
routing model' met these criteria and also those of hi^ remote sen- 
sing input potential and ocmputational sinp>licity. 
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5.3) EE^drologic Analysis of LSNDSftT Imagery 

The final task was cdmed at using remote sensing directly to deter- 
mine hydrologic information content of the LANDSAT bands and to at- 
atn^Jt to sir^y extract the necessary hydrologic data. Initially, 
the Oxford and Blacksburg test basins were examined. It was found 
that information sufficient to (fetermine several of the important 
inputs to the model could be determii^d from LANDSST data using re- 
latively uncoT^jlicated visual techniques. Jtoteover, it was deter- 
mined that single bands contain useful but different data and hypo- 
thesised that a ocmposite image \iDuld optimize the information value. 
This l^pothesis was confimed through in-(fepth analysis of the Middy 
Branch basin. Surface features of the watershed were identified, 
measured, and checked against aerial photogr^ahs ground truth. The 
results showed inverttory errors to be well within acceptable limits 
for modelling and useful for direct ccmtputation of model parameters. 


APPE3SDIX A 


STAnSUCftL AEggjYSIS OF FOfECAST ERIPR Bffi>K3VEMEWP 


A set of tests was run to determine if the observed ijnpravanstits were 
statistically significant. The proce(3ure to be used is to first es- 
tablish the equality or inequality of variances among the four models 
and given these results, to determine the significance of the differ- 
ence among the means. First, the Oochran test for equality of vai- 
anoe was performed to test the nvill hypotl^sis that 

2 ^2 

®peak rate “ °^SCS, Cook, Eational 
, against the alternative that the variances are unequal. 


The e^ropriate statistic is as follows! 

2 


G = 


n ; 

I 


i«l 




Where: G = test statistic number 

= sanple variance of sanple 

k = numter of treatments {equal to 2 in this case 
since the models were tested in pairs) 


The process for ccwparing the variances of the peak rate and S.C.S. 
models is given below. 



peak irate 


6.691.24 

1.332.25 


r - 6,691.24 
^ "■ 8,023.45 


0.83 


The table value of G at the 95% cxanfidence level is 0.7657. 
test table 

The null hypothesis for the peak rate and SCS formulas can therefore be 
rejected, thus accepting the alternative hypothesis. 


a 2 

“scs > a , 

peak rate 

Similarly, the equality of variance of the peak rate and Cook's model was 


tested. 


peak rate 


Cooks 


1,332.25 

22,290.49 


22,290.49 

23,622.74 


0.94 


^test ^ ^table 


( - 0,7f-57 at 95% confidence level) 


we reject the null hypothesis, and oonclude that 


CX** ^ ZL 

peak rate “Oboks 

Further, the inequality of variances between peak rate and Rational was 
established 


peak rate 
S*' 

Rational 


1,332.25 

26,634,24 


o 


G 


26.634,24 _ 

27,966.49 


o o 


*^test ^ Sable rejc:cted tha null fe/pcDthesis in favor of the alter- 
native hypothesis, and conclude that 

2 ^2 
peak rate ^ °T>ational 

All the above analyses enable us to conclude that tha variance of the 
peak rate iroctel is less than tiie vari£uic''s of the other models. The 
reduction of variance observed is therefore statistically significant 
at the 95% confidence le‘/ei. 


Knowing that the variances between the peak rate model and the other 
three are not equal defines the appropriate test for equality of the 
means. In this series, the peak rate model was tested indi\rLdually 
against the three other models for significance of difference between 
tha means. The test statistic losed was: 


T = 


S 2 ^2 

1 , "2 
4 * 




n. 


T = test stotistic 

^ 2 ~ errors of the two test samples 

2 2 

^1 ' ^2 ~ variances of the sanpl er 

Hi, ^2 ~ of each sanple 


In the first the peak rate model was crrpared with the SCS for- 


peak rate 


= 54,1 


“PCS 


- 62.5 


mula: 


6,691.24 


^ « 1,332.25 
peak rate ' 


^ peak rate ~ 


T = 


SCS 

62.5 - 54.1 


« 31 


6,691.24 , 1,332.23 
31 


D.s; 


The OTmber of degrees of ^Ireedom were ocaaputed by the formula: 



41.49 


« 42 degrees of freedom 

The value of the table statistic at the 90% confidence le\iel for 42 
degrees of freedon is less than 1.311 


o fT> < 4 - 

• * test ralde 


The hypothesis of equality of means must therefore be accepted showing 


the peak rate model, errors to be ecfual to the SCS nv:x3e] . 


The second test ccatpared the peak rate and Cook irodels. 


^Ooo]c 

.2 

* Oook 
*^0ook 


99.2 

22,290.49 

31 
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tha same equation as previ.ously given: 

T « 1.63 
V 34 

The table value of at 34 degx-ees of freedom is less than 1.311 

test ' 0.05 o‘o Peject tl'^s null hypothesis and con- 

clude that the mean pealt rate model anx r is significantly less than 
that of the Oook formula. 


Finally, the peak rate model was (;;-3mpar.ed with the national fomwla. 


Therefore, 


Rational 

I^ational 


86.7 

26,634,24 


^^tional 


T 1.09 
V - 33 


The significance of error- impxovenent oould not be established at thi 


confidence level. 
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